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We study the growth of perturbations in a uniformly collapsing cloud of self-gravitating Brownian 
particles. This problem shares analogies with the formation of large-scale structures in a universe 
experiencing a "big-crunch" or with the formation of stars in a molecular cloud experiencing grav- 
itational collapse. Starting from the barotropic Smoluchowski-Poisson system, we derive a new 
equation describing the evolution of the density contrast in the comoving (collapsing) frame. This 
equation can serve as a prototype to study the process of self-organization in complex media with 
structureless initial conditions. We solve this equation analytically in the linear regime and com- 
pare the results with those obtained by using the "Jeans swindle" in a static medium. The stability 
criteria, as well as the laws for the time evolution of the perturbations, are different. The Jeans 
criterion is expressed in terms of a critical wavelength Xj while our criterion is expressed in terms 
of a critical polytropic index 74/3. In a static background, the system is stable for A < Aj and un- 
stable for A > Aj. In a collapsing cloud, the system is stable for 7 > 74/3 and unstable for 7 < 74/3. 
If 7 = 74/3, it is stable for A < Aj and unstable for A > Xj. We also study the fragmentation 
process in the nonlinear regime. We determine the growth of the skewness, the long-wavelength 
tail of the power spectrum and find a self-similar solution to the nonlinear equations valid for large 
times. Finally, we consider dissipative self-gravitating Bose-Einstein condensates with short-range 
interactions and show that, in a strong friction limit, the dissipative Gross-Pitaevskii-Poisson system 
is equivalent to the quantum barotropic Smoluchowski-Poisson system. This yields a new type of 
nonlinear mean field Fokker-Planck equations including quantum effects. 



I. INTRODUCTION 

This paper continues our theoretical investigations (see 
[l[ for a short review) of the barotropic Smoluchowski- 
Poisson (BSP) system ([T|)-(I21). For an isothermal equa- 
tion of state p = pkBT/m, these equations describe 
the mean field dynamics of a dissipative gas of self- 
gravitating Brownian particles in the overdamped limit 
^ — ?> -l-oo where inertial effects are neglected. For the sake 
of generality, we shall consider an arbitrary barotropic 
equation of state p = p{p) so as to take into account 
anomalous diffusion or short-range interactions. This 
type of equations appears in different contexts (with dif- 
ferent interpretations) such as planetary formation in 
the solar nebula 0, chemotaxis of biological popula- 
tions 3, colloids at a fluid interface driven by attractive 
capillary interactions [1] and nanoscience d, Q. They 
also provide a simplified model of gravitational dynam- 
ics which displays many features common with ordinary 
self-gravitating systems such as collapse and evaporation 
[H . An interest of this model is that many analytical re- 
sults can be derived since the inertia of the particles is 
neglected. For these reasons, the theoretical study of the 
BSP system is of considerable importance. 

When studying equations of that type, it is natural 
to investigate first the dynamical stability of a spatially 
homogeneous distribution. However, in the gravitational 
case, we are rapidly confronted to the well-known prob- 
lem that a spatially homogeneous distribution is not a 
steady state of these equations. To circumvent this dif- 
ficulty, one possibility is to advocate the "Jeans swin- 
dle" as in similar problems of astrophysics based on the 
barotropic Euler-Poisson system [S-Q- This procedure 



allows one to perform the stability analysis in a simple 
manner. Although this first-step study certainly sheds 
light on the problem, it is however not satisfactory and 
must be improved. Another possibility is to study the 
dynamical stability of spatially inhomogeneous steady 
states. This has been done in [HI for isothermal 
spheres, in [l^, [l!] for polytropic spheres and in for 
an arbitrary equation of state by extending the methods 
developed in astrophysics |15l - [3(]| |. Finally, a third pos- 
sibility is to construct time-dependent solutions of the 
BSP system. Solutions describing the collapse, the post- 
collapse and the evaporation of the system have been con- 
structed in [l0l - [l^ . l3TI - [37j by looking for self-similar, or 
close to self-similar, solutions. A fully analytical solution 
of the BSP system has been obtained in [3^ in the cold 
case where the pressure can be neglected. [Similar studies 
have been previously done in as trophy sics based on the 
barotropic Euler-Poisson system |39l - [45j . The results are, 
however, different because the barotropic Euler-Poisson 
system takes into account the inertia of the particles but 
neglects dissipation while the BSP system is obtained in 
a strong friction limit in which the inertia of the particles 
can be neglected. The barotropic Euler-Poisson system 
and the BSP system correspond therefore to opposite lim- 
its (weak and strong dissipation).] 

Among these solutions, there is a very simple one, 
namely the collapse of a homogeneous sphere. It is found 
[sst that the sphere undergoes a finite time singularity, 
resulting in a vanishing radius and an infinite density 
in a finite time. The question that naturally emerges 
is whether this time-dependent solution is stable with 
respect to small perturbations and, if not, how the per- 
turbations grow. This amounts to considering the de- 
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velopment of density fluctuations in a time-dependent 
collapsing cloud. In a sense, this is the right formulation 
of the Jeans problem without "swindle" . We may note 
the analogy with the formation of stars in astrophysics 
through the collapse and the fragmentation of a molec- 
ular cloud [H, 13 53] • While the cloud collapses under 
its own gravity, it also develops numbers of irregulari- 
ties that grow and finally form stars. We may also note 
the analogy with the formation of large-scale structures 
in cosmology HI] except that the universe is expanding 
[s?! while, in our problem, the cloud is contracting. In 
a sense, this corresponds to a universe experiencing a 
"big-crunch" . To perform our study we shall use many 
methods introduced in cosmology [il,!!^]- There wih be, 
however, two main differences: (i) our system is contract- 
ing instead of expanding; (ii) inertial effects are neglected 
in our model. The second point is a huge simplification 
with respect to cosmology that allows us to solve the lin- 
ear stability problem fully analytically. We can therefore 
easily follow the growth of perturbations in the linear 
regime for an arbitrary equation of state. Furthermore, 
we derive an exact system of nonlinear equations [see Eqs. 
([M|) - (|35p ] describing the formation of clusters on longer 
times. These equations are simpler than those considered 
in cosmology and could serve as a prototype to study the 
process of self-organization in complex media with struc- 
tureless initial conditions. 

The paper is organized as follows. In the first part of 
the paper (Sees. IIIllIV|) . we consider the classical BSP 
system (IT|)-(I21). In Sec. |lll we study the dynamical sta- 
bility of an infinite homogeneous distribution in a static 
frame by invoking the "Jeans swindle" . We find that the 
system is stable if and only if (iff) the wavelength of the 
perturbation is smaller than the Jeans length Aj. The 
perturbations decay exponentially rapidly if A < A,/ and 
grow exponentially rapidly if A > A,/. In Sec. IIIIl we 
derive a new set of equations [see Eqs. ([M]). ([55]) and 
(P^ j describing the evolution of the density contrast in 
a homogeneous collapsing cloud. We study the develop- 
ment of perturbations in the linear regime and obtain an 
analytical solution for the time evolution of the Fourier 
modes of the density contrast. For a polytropic equa- 
tion of state, we find that the evolution depends on the 
value of the polytropic index 7. There exists a criti- 
cal index 74/3 = 2(d — l)/d, already encountered in our 
previous studies jl2Hl4l [ and well-known in astrophysics 
[1^ [2l|, such that the system is stable for 7 > 74/3 and 
unstable for 7 < 74/3. The effect of the Jeans length 
Xj{t), which is a function of time, manifests itself only in 
the initial stage and disappears as time goes on. In the 
unstable case, the perturbation is always growing when 
A > Aj(0) while it starts decaying, before finally growing, 
when A < Aj(0). We find that the perturbation grows 
algebraically for large times. This is in sharp contrast 
with the exponential growth predicted by the "naive"' 
Jeans analysis in a static medium. In the stable case, 
the perturbation is always decaying if A < A j (0) while it 
starts growing, before finally decaying, when A > A,/(0). 



For 7 = 74/3, the Jeans length is constant. The pertur- 
bation grows for A > Aj and decays for A < Aj. This is 
similar to the result of the Jeans analysis except that the 
evolution is algebraic instead of exponential. In Sec. IIVI 
we study the nonlinear regime and derive several analyti- 
cal results. In particular, we determine the growth of the 
skewness and the long- wavelength tail of the power spec- 
trum. We also find an exact self-similar solution to the 
nonlinear equations valid for large times. In the second 
part of the paper (Sec. |V]), we include quantum effects 
in our model. Specifically, we consider a self-gravitating 
Bose-Einstein condensate (EEC) in the presence of dissi- 
pative effects. Self-gravitating BECs have been proposed 
recently as a model of dark matter in cosmology [sif . 
Starting from a dissipative form of the Gross-Pitaevskii- 
Poisson system and using the Madelung transformation, 
we derive the quantum damped barotropic Euler-Poisson 
system. Considering a strong friction limit, we obtain 
the quantum barotropic Smoluchowski-Poisson system 
(|120p - (|12ip . This corresponds to the classical barotropic 
Smoluchowski-Poisson system (H])-© with an additional 
term called the Bohm quantum potential. We generalize 
our previous analysis in this more general context. 

II. SMOLUCHOWSKI-POISSON SYSTEM IN A 
STATIC FRAME 

A. Barotropic Smoluchowski-Poisson system 

We consider the barotropic Smoluchowski-Poisson 
(BSP) system 

e|^ = V.(Vp + pV$), (1) 

A$ = SdGp, (2) 

where p{r, t) is the spatial density of the particles, $(r, t) 
is the gravitational potential, p(r, t) is the pressure, ^ 
is the friction coefficient and Sd is the surface of a unit 
sphere in d dimensions. Basically, this equation describes 
the competition between pressure effects and gravity. To 
close the system, we assume that the pressure is related to 
the density by an arbitrary barotropic equation of state 
p = p{p). The case of (ordinary) Brownian particles cor- 
responds to an isothermal equation of state 



where the constant T is the temperature of the bath and 
TO the mass of the particles (fc^ is the Boltzmann con- 
stant). Another very useful, and popular, equation of 
state is the polytropic equation of state 

p = Kp\ 7=1 + -, (4) 
n 

where K is the polytropic constant and 7, or n, is the 
polytropic index. The isothermal case is recovered for 
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7 = 1, i.e. n — >■ cxD, and K = kBT/m. The velocity 
of sound Cs is defined by the relation = p'{p)- For an 
isothermal equation of state = ksT/m is constant and 
for a polytropic equation of state = K^p'~'^^ depends 
on the density. 

The Smoluchowski-Poisson system with an isothermal 
equation of state ([3]), corresponding to self-gravitating 
Brownian particles, has been derived in [52| . It is valid 
in a mean field approximation that becomes exact in a 
proper thermodynamic limit N — ^ +00 such that the 
normalized temperature rj = /SGMm/ R'^~'^ is unity {R 
denotes the size of the system and (3 = l/ksT the inverse 
temperature). The barotropic Smolchowski-Poisson sys- 
tem with an arbitrary equation of state p{p) has been 
derived in [53l [H^l from a notion of generalized thermody- 
namics and in ^531 from the dynamic density functional 
theory (DDFT). In the two approaches, the nonlinear 
pressure arises from short-range interactions between the 
particles but the physical arguments leading to Eq. ([T]) 
are different. 



B. Jeans- type instability 

A steady state of the BSP system ([I])-© satisfies the 
relation 



Vp + pV$ = 0, 



(5) 



which can be viewed as a condition of hydrostatic bal- 
ance. Since p = p{p), this relation can be integrated 
to yield p = p{^)- Using the Poisson equation ([5]), we 
obtain the differential equation 
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which has to be solved with appropriate boundary con- 
ditions. We note that self-gravitating Brownian parti- 
cles described by the barotropic Smoluchowski-Poisson 
system have the same equilibrium states as barotropic 
stars described by the barotropic Euler-Poisson system 
[^,[13. However, their dynamics is different since the mo- 
tion of Brownian particles is overdamped. The barotropic 
Smoluchowski equation ([T]) can be written 



,dp ^ f ^SF' 



(7) 



where F is the free energy functional 

F = ^Jp'i>dr + JpJ'^^dpidr. (8) 
The BSP system satisfies the i/-theorem 

Therefore, the BSP system relaxes towards the state that 
minimizes the free energy at fixed mass. This corre- 
sponds to the condition of hydrostatic equilibrium 



In this sense, the BSP system is consistent with the sta- 
tistical equilibrium state in the canonical ensemble. By 
contrast, the functional dH) is conserved by the barotropic 
Euler-Poisson system so that there is no relaxation to- 
wards equilibrium in that case 

A natural problem concerns the dynamical stability of 
a steady state of the BSP system. We shall restrict our- 
selves to linear dynamical stability. To that purpose, we 
consider a small perturbation Spirit) around a steady 
state p{r) of Eqs. ([I])-© and study its dynamical evo- 
lution. In the general case, the linearized BSP system 
is 



dSp 
'dt 



V ■[V{cl{p)5p) + 8pV^ + pV5<^], (10) 



= SdGSp. 



(11) 



The dynamical stability of box-confined isothermal and 
polytropic configurations [s^ has been studied in [lol - 
[l3| . For example, for the isothermal equation of state 
in d = 3, it is found that the system is stable if 
rj = jSGMm/R < t]c ~ 2.52 and if the density con- 
trast n EE p{0)/p{R) < Uc 32.1; it is unstable 
if one of these two conditions is not met [2^. Simi- 
larly, for a polytropic equation of state in d = 3, it is 
found that incomplete polytropes with 7 < 6/5 (n > 5) 
are stable iff their inverse polytropic temperature r/ = 
M[47rG/if(n-Hl)]"/("-i)/47ri?("-3)/("-i) and their den- 
sity contrast TZ = p{0)/p{R) are sufficiently small (the 
exact values of r]c{n) and TZc(n) depend on the index 
n) [2^. The dynamical stability of a general equation 
of state p(p) for self-confined configurations has been 
studied in [ij] by adapting the methods of astrophysics 
[isl [2]| . For a polytropic equation of state in d dimen- 
sions, it is found that a spatially inhomogeneous self- 
confined configuration is stable iff 7 > 74/3 = 2{d— l)/d. 
When the steady state is unstable, it usually collapses 
or evaporates depending whether diffusion (pressure) or 
attraction (gravity) prevails [l[. 

Let us now consider the case of a spatially homoge- 
neous distribution. In that case, Eqs. pO)) - ((TT|) reduce 
to 



d6p 
^ dt 



lASp + SdGpSp, 



(12) 



where cj. = c^ip) is constant. Considering for simplicity 
an infinite system [89| and decomposing the perturba- 
tions in plane waves of the form g*(k r-Ljt)^ obtain the 
dispersion relation 



SdGp. 



(13) 



This relation, like the original Smoluchowski equation 
([1]) . clearly shows the competition between the attractive 
gravitational force and the repulsive pressure (assuming 



> 0). There exists a critical wavenumber 



SdGp 



1/2 



(14) 
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corresponding to the Jeans wavenumber in astrophysics 
However, contrary to the astrophysical problem 
based on the barotropic Euler-Poisson system, the evo- 
lution of the perturbations is different in our case. This 
is because the Smoluchowski-Poisson system is an over- 
damped model where inertial effects are neglected. Here, 
the perturbations evolve like e'''* with an exponential rate 
7 = {SdGp — c^fc^)/^. For k < kj, the system is un- 
stable and the perturbations grow with a growth rate 
7 > 0; for fc > fcj the system is stable and the pertur- 
bations decay with a damping rate 7 < 0. Note that the 
growth rate is maximum for k = (corresponding to an 
infinite wavelength A = 27r//s — > +00) and its value is 
"fmax = SdGp/^. In conclusion, Jeans' stability criterion 
is k > kj (i.e. X < Xj ). 

Although this stability analysis, like the original Jeans 
study in astrophysics 0i Hi is simple and sheds light 
on the problem, it however encounters a serious prob- 
lem. Indeed, if p is uniform (in a finite or infinite do- 
main) , the pressure gradient in Eq. ([5]) vanishes implying 
that the gravitational field must also vanish (V<f> = 0). 
Then, the Poisson equation © cannot be satisfied for 
p > 0. This implies that there is no steady state of 
Eqs. (O-® with uniform density. Therefore, when us- 
ing the perturbation equation (jl2p . we make what is tra- 
ditionally called the "Jeans swindle" @. There are dif- 
ferent ways to circumvent this difficulty, (i) The first is 
to restrict ourselves to spatially inhomogeneous steady 
states that are solution of Eqs. ([S])-®. (ii) The sec- 
ond is to consider a uniformly rotating system. In that 
case, if we introduce an effective gravitational poten- 
tial $6// = — (n X r)^/2 in the rotating frame, the 
Poisson equation is replaced by A$e// = 47rGp — 2^2^ 
and there exists spatially homogeneous distributions with 
p = ^ 12-kG. (iii) The third possibility is to modify the 
Poisson equation. This is justified for some physical sys- 
tems. For example, in the context of chemotaxis 0, the 
Poisson equation is replaced by a generalized field equa- 
tion of the form x9f $ = A*!" — k\^ — SdGp where kj^ is 
a screening length. This leads to the general Keller-Segel 
model. In that case, a spatially homogeneous distribution 
is a steady state of the equations with $ — —SdGp/k\. 
In some approximations, the generalized field equation 
= A<I> — A:|j4> — SdGp is replaced by the screened 
Poisson equation A4> — A:|,<I' = SdGp or by the modified 
Poisson equation A<I> = SdG{p—'p) where p is the average 
density. Again, a spatially homogeneous distribution sat- 
isfying $ = —SdGp/kj^ or p = p respectively is a steady 
state of these equations. Therefore, there is no "Jeans 
swindle" in the chcmotactic problem [56l - l58| . In grav- 
ity models, we could introduce a cosmological constant 
A and replace the Poisson equation by A$ = SdGp — A. 
In that case, the particular homogeneous density distri- 
bution p = A/SdG is a steady state of the equations fooj . 
However, this justification does not apply to a uniform 
mass of gas of different density. Furthermore, this homo- 
geneous distribution is unstable as shown in Appendix 
IB] (iv) Kiessling [6l| has provided a vindication of the 



Jeans swindle. He argues that, when considering an infi- 
nite and homogeneous distribution of matter, the Poisson 
equation must be modified so as to correctly define the 
gravitational force. He proposes to use a regularization 
of the form A$ — kj^^ = AirGp where fcjj is an inverse 
screening length that ultimately tends to zero {ku — > 0), 
or a regularization of the form A$ = AttG{p — p) where 
p is the mean density. In his point of view, the modifi- 
cation of the Poisson equation is not a swindle but just 
the right way to make the problem mathematically rigor- 
ous and correctly define the gravitational force. However, 
this argument applies only to an infinite system. For a 
finite system, a spatially homogeneous distribution is not 
a steady state of the BSP equations. 

Remark: a form of Jeans criterion can be recov- 
ered for box-confined configurations. For example, for 
finite isothermal spheres in d ~ 3, introducing the 
mean density p = 3M/47ri?'^, the stability criterion 
■q = (3GMm/R < rjc ~ 2.52 ^ can be rewritten 
R < Rj = (377c/47rG'pm/3)^/^ where Rj is similar to 
the Jeans length Xj = 2'K{kBT / AttG pmY^^ correspond- 
ing to Eq. p^ . Similarly, for finite polytropic spheres 
with n > 5, the stability criterion 77 = MlAirG / K{n + 
lyn/in-i) i^^p,(n-i)/{n-i) ^ ^ ^.^^ rewrit- 

ten R < Rj = [K{n+ l)0(n)/47rGnp("-i)/"]i/2 with 
(?!)(n) = n[377c(n)]("^^'/" where Rj is similar to the Jeans 
length Xj = 27r[X(n -f l)/47rG7ip("~i)/"]i/2 correspond- 
ing to Eq. (fT4| . 



C. Homogeneous collapsing cloud 

It is clear that a finite spatially homogeneous solution 
of the BSP system will collapse under its own gravity. 
Similarly, a finite spatially homogeneous solution of the 
Eulcr-Poisson system is not steady and will evolve in 
time. However, since the Euler-Poisson system is an in- 
ertial model, depending on its energy, the system will 
either expand (positive energy) or collapse (negative en- 
ergy) . This is at the basis of the Newtonian cosmological 
models of McCrea & Milne [gJI where the universe is 
viewed as a homogeneous self-gravitating sphere of ra- 
dius a{t) [9H. In cosmology, we can have open (always 
expanding) or closed (expanding then contracting) mod- 
els of the universe with an intermediate ca se p rovided 
by the Einstcin-de Sitter (EdS) universe mm- In our 
problem, since the velocity of the particles is directly pro- 
portional to the gravitational force (overdamped limit), a 
homogeneous cloud can only collapse under its own grav- 
ity. In the analogy with cosmology, this corresponds to 
a "big-crunch" . This is also similar to the gravitational 
collapse of a molecular cloud in relation to the process of 
stars formation [3^,|4l|43. 

It is easy to construct a time-dependent spatially ho- 
mogeneous solution of the BSP system. For a spatially 
homogeneous distribution of the form 

p(r,i)=pfc(t), (15) 
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the BSP system reduces to 

dpb SdG 2 

The pressure p does not appear in the theory of the ho- 
mogeneous model since it enters Eq. ([1]) only through 
its gradient. Therefore, we recover the equation of a cold 
system (p = 0) considered in [sl]- If we consider the 
spherically symmetric collapse of a uniform sphere of ra- 
dius a(i), the conservation of mass imposes 



a 



(17) 



Solving Eq. ([T5|) and using Eq. ((T7)) . we find that 
the density and the radius of a collapsing homogeneous 
sphere evolve according to 

PbiO) 



Pb{t) = 



ait) = a(0) 1 



l/d 



where 



SdG PbiO) dGM ■ 



(18) 



(19) 



(20) 



This solution leads to a finite time collapse. At t = , all 
the mass is in a Dirac peak at r = (big crunch). This 
solution is the counterpart of the homogeneous solution 
found by Hunter [si], Mestel ^ and Penston [U in the 
inertial case. 

We can also obtain these equations from the La- 
grangian motion of a fiuid particle as in jssj by using the 
fact that the Smoluchowski equation ([T]) with Vp ~ 
can be interpreted as a continuity equation with a veloc- 
ity field 

:V$. (21) 



dt 



Applying this equation at the border of the sphere, and 
using the Gauss theorem, we obtain 



da 

ctt ' a- 



GM 



-SdGpba, 
d 



(22) 



which after integration returns Eq. (|T9)) . In the analogy 
with cosmology, Eqs. (fTT]) and ([22]) are the counterparts 
of the Friedmann equations when pressure effects are ne- 
glected. We note, however, that the Friedmann equations 
arc second order in time while Eq. (|22p is first order. 



III. 



SMOLUCHOWSKI-POISSON SYSTEM IN A 
COLLAPSING FRAME 



This process of fragmentation, driven by gravity, ulti- 
mately leads to the formation of localized dense clusters 
similar to galaxies in the universe or stars in a molecular 
cloud. Like in cosmology, we shall work in the comoving 
frame (here the collapsing frame). To that purpose, we 
set 



a(t)x. 



(23) 



where ait) is the scale factor (in our problem it corre- 
sponds to the radius of the sphere). Let us first write 
the Poisson equation ([2]) in the collapsing frame. For the 
background distribution (|15p . the Gauss theorem reads 



d$b GMir) 1 a 
-J- = w- 1 = -iSdGpbr = -C-r, (24) 
ar r°- ^ a a 

where A/(r) is the mass enclosed within the sphere of ra- 
dius r and we have used Eq. ((22|) to obtain the last equal- 
ity. The background gravitational potential is therefore 



^ 1^09 1 2 
$fc = -T^S.-r- = --£,aax\ 
2 a 2 

In terms of the new potential = •I' — $b, i.e. 

1 „ 
</) = $-!- -^aax , 

the Poisson equation ([2|) becomes 

A^^SdGa^ip^Pb). 



(25) 



(26) 



(27) 



In this expression, the Laplacian is defined with respect 
to the variable x and we have used Eq. ([22|) to obtain 
the second term on the r.h.s. 

We now write the barotropic Smoluchowski equation 
in the collapsing frame. We first note that 



I) p(rM.),.) = |-^x.Vp. (28) 



Therefore, the barotropic Smoluchowski equation ([1} be- 
comes 

^--x.Vp=-^V-(Vp + pV<I>), (29) 
ot a t,a'^ 

where the derivatives are taken with respect to x. Intro- 
ducing the new potential ([26)) . we obtain after simplifi- 
cation 

^+d-p^^V ■iVp + pV<t>). (30) 
ot a 

Like in cosmology, it is convenient to write the density 
in the form 



A. The fundamental equations 

We shall now study the instability of the collapsing ho- 
mogeneous sphere and the development of irregularities. 



p = Pbit) [l + 5(x,t)] 



(31) 



where pb oc l/a"^ according to Eq. pT)) and (5(x,i) is the 
density contrast Substituting Eq. in Eq. ([50]) . 
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we obtain the barotropic Smoluchowski equation in the 
collapsing frame 



Pb 



as 



1 



dt 



2V-(Vp + /jfc(l + ,5)V0). 



Since p = p{p), this can be written 



dS 



1 



- = -V.(c^V.+ (l + .)V0), 



(32) 



(33) 



where = p'{p) = p'{Pb{^ + ^)) is the square of the 
velocity of sound in the evolving system. It generically 
depends on position and time (it is constant only for an 
isothermal equation of state). Finally, substituting Eq. 
(|3T1) in the Poisson equation ((27)) . we obtain our final 
system of equations 



^dt a? 



Act) = SdGpbSa^, 



(34) 



(35) 



where pb{t) and a{t) are given by Eqs. ([T5)) and p^ . If 
we measure the evolution in terms of a instead of t, we 
can rewrite Eq. p4p in the form 



da 



SdGpbW 



rV-(c2V(5+(l + <5)V</.), 



(36) 



where we have used Eq. (I22|) . For an isothermal equation 
of state, Eq. ([34]) becomes 



Ot a'^ \ m J 

(37) 

where we have used the Poisson equation (|35l) to get the 
second term on the r.li.s. For a cold system = 0, we 
get 



= SaGptS{l 

at 



S) + \vS ■ V0. (38) 



Equations p4p -([35 |) could form a prototypical model to 
study the formation and the growth of structures in com- 
plex media with structureless initial condition. As we 
shall see, they lead to a process of fragmentation and 
self-organization. Furthermore, they are simpler to study 
than the equations used in cosmology since there 
is no inertia in our model. Indeed, Eqs. (p4|) -([35 |) are 
parabolic while the equations used in cosmology are hy- 
perbolic. We therefore expect to obtain simpler, hope- 
fully analytical, results like in the study of the BSP sys- 
tem (P-dl]). 



flow ([T5|) . We shall now study the dynamical stability of 
this solution. If we consider small perturbations 5 <C 1, 
<C 1 around that state, we obtain the linearized equa- 
tion 



85 
dt 



;AS + SdGpbS, 



(39) 



where ~ p'ipt) is the square of the velocity of sound in 
the background flow. It depends only on time. Equation 
(1391) is the counterpart of the equation derived by Bonnor 
[48| in cosmology. Expanding the solutions of Eq. ([M)) 
in Fourier modes of the form S{x, t) = 5k(i)e*''", we find 
that 



^5 



- SdGpb U = 



(40) 



where, for convenience, we have noted 6 for S]^. From 
this equation, we can define a time-dependent Jeans 
wavenumber in the comoving frame 



_ / SdGpba^ 



1/2 



(41) 



Note that the proper wavenumber obtained by writing 



(5 oc e 



is k*i 



{SdGpb/c 



,1/2 



The evolution of 



kj{t) depends on the equation of state. For an isothermal 
equation of state for which = ksT/m, we can write 
kj = Kja-(''-2)/2 where kj = {SdGpba'^m/kBTy/^ is a 
constant according to Eq. (fT7| . The Jeans length Aj = 
2TT/kj behaves like a^'^'^'^/^. For d > 2 (resp. d < 2), it 
decreases (resp. increases) with time so that the system 
becomes unstable at smaller and smaller (resp. larger 
and larger) scales as the cloud collapses. For d = 2, the 
Jeans length is constant: kj = kj. For a polytropic equa- 
tion of state for which = KjpJ~^ , we can write kj = 

^^^(d7-2(d-l))/2 ^ j ^ [5dG(p6a'*)2-T/if7]l/2 ig 

constant according to Eq. (fTT]) . The Jeans length be- 
haves like Aj oc a(2(rf-i)-'i7)/2, For7 < 74/3 = 2(d-l)/d, 
the Jeans length decreases with time and, for 7 > 74/3, 
it increases with time. For 7 = 74/3, the Jeans length is 
constant: kj = kj. The evolution of the proper Jeans 
length A J is different but the evolution of the comoving 
Jeans length Aj is better suited to our problem. 

Remark: our analysis exhibits a critical index 74/3 = 
2{d — l)/d, the same as the one determining the dy- 
namical stability of spatiall y in homogeneous polytropic 
spheres p"3 - [l^ (see also [H,[2l| in astrophysics) . It is in- 
teresting to note that this critical index does not appear 
when the Jeans problem is (incorrectly) formulated in a 
static homogeneous background (see Sec. IIIB[) [o^ while 
it appears when the Jeans problem is properly formu- 
lated in a collapsing background. This will have impor- 
tant consequences in the stability problem (sec below). 



B. Linearized equations 



C. Solution of the equation for the density contrast 



Of course, 5 = </> = is a particular solution of 
Eqs. ([M)) - (|55|) corresponding to the pure background 



The evolution of the perturbations is more complicated 
to analyze in a collapsing homogeneous medium than in 
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a static homogeneous medium (and turns out to be very 
different). An interest of the present model is that the 
equation (|40|) for the density contrast in the hnear regime 
can be solved analytically for any equation of state while 
the corresponding equation in cosmology can be solved 
analytically only in particular cases (this is because it is a 
second order differential equation in time while Eq. (|40p 
is a first order differential equation) |i49j. Measuring the 
evolution in terms of a instead of t and using Eq. (|22]) . 
we can rewrite Eq. (j40p in the form 



da a \SdGpba? 
Its general solution is 



(42) 



(43) 



For a cold gas = in which the pressure is zero, Eq. 
dH]) reduces to 



d5 d 



5 = 0, 



da a 
and its solution is 



(44) 



(45) 



In that case, the homogeneous cloud is unstable at all 
scales and we can write in full generality (5(x,a) = 
D{x)/a'^ or S{x,t) = (5(x, 0)(l-i/**)"^- We note that the 
growth of the perturbation is algebraic while the static 
Jeans study predicts an exponential growth with a rate 
7 = SdGp/£, (see Sec. Ill Bp . This is a very important dif- 
ference. Assuming now a polytropic equation of state and 
introducing the notations of Sec. IIII Bl we can rewrite 
Eq. (|42]) in the form 



d6 



da a V K^a'i-i-^i'^-^') 



- 1 U = 0. 



(46) 



If 7 7^ 74/3 (i-C. n ^ ^-s), the solution of this equation 



.2 2{d-l)-d^ 




FIG. 1: Evolution of the perturbation S{a) for 7 > 4/3 {d = 3, 
ao = 1 and So = 1). We have taken 7 = 2 which corre- 
sponds to a BEG with quartic self-interactions in the TF ap- 
proximation (see Sec. IV C|) . The system is asymptotically 
stable. For k < fcj(O), the perturbation grows before de- 
caying. The maximum occurs at a = {k/kj)^^^^''~^^ and its 
value Smax oc fc~'^/(-3T~'*) is maximum for small k (large wave- 
lengths). 
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1,1," 



(47) 



1/a 



FIG. 2: Evolution of the perturbation 5(a) for 7 < 4/3 
(d = 3, ao = 1 and (5o = 1). We have taken 7=1 
which corresponds to a classical self-gravitating Brownian gas 
(isothermal distribution). The system is asymptotically un- 
stable. For k > fcj(O), the perturbation decays before grow- 
ing. The minimum occurs at a = (fc/fcj)~^/'*~'^^' and its 
value Smin. OC is minimum for large k (small wave- 

lengths). 



This includes the isothermal equation of state (7=1) 
with d 7^ 2 as a particular case 



1 



8(a) oc e'''-^'"/ 



(48) 



For short times, writing a = ao(l — e) with e ^ 1, we 
obtain 



(5(ao) 



1 + l 



fcj(O)^ 



(49) 



This equation clearly shows the initial effect of the Jeans 
wavelength. Indeed, the perturbation starts to grow 



when k < kj{Q) and starts to decay when k > fc,/(0). 
This result, which is independent on the value of 7, is 
similar to the one obtained from the static Jeans study 
(see Sec. Ill Bp . Indeed, for short times, the system docs 
not "see" that the background is contracting. By con- 
trast, for large times, the effect of the Jeans scale disap- 
pears and the evolution of the system is only controlled 
by the value of 7. For 7 > 74/3, S{a) — >■ very rapidly 
as a ^ 0. The system becomes asymptotically stable to 
all wavelengths. For 7 < 74/3, S{a) cx a~'^ — > +00 as 
a — >■ 0. The system becomes asymptotically unstable to 
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FIG. 3: Evolution of the perturbation S{a) for 7 = 4/3. The 
system is stable for k > kj and unstable for k < kj. We have 
taken ao = 5o = 1- 



cases, an imaginary velocity of sound arises from a neg- 
ative pressure due to attractive short-range interactions 
(see Sec. IV Cp . In that case, we find that the perturba- 
tions S{a) always grow. Therefore, a negative value of 
can enhance the gravitational instability and make the 
system very unstable. 

IV. NONLINEAR GRAVITATIONAL 
CLUSTERING 

In this section, using an analogy with cosmology, we 
provide analytical results valid in the nonlinear regime. 

A. Second order perturbation theory for Sp/pi,: 
skevirness 



all wavelengths and behaves like in a cold gas. There- 
fore, the system is asymptotically stable for 7 > 74/3 and 
asymptotically unstable for 7 < 74/3. This stability cri- 
terion is very different from the Jeans stability criterion 
in a static background (sec Sec. Ill B[) where the value 
of the index 7 docs not play any crucial role. It is, on 
the other hand, compatible with the stability criterion of 
self-confined inhomogeneous polytropic spheres [l^ . 
If 7 = 74/3 (i.e. n = ns), the solution of Eq. (^5)) is 



(50) 



This includes the isothermal equation of state (7=1) 
with d = 2 as a particular case. We see that the Jeans 
scale kj ~ kj now plays a crucial role at all times. For 
k < Kj, the perturbation increases in time and the sys- 
tem is unstable (for k <^ kj, the perturbation grows like 
in a cold gas). For k > kj, the perturbation decreases 
in time and the system is stable. These results are sim- 
ilar to those obtained in the static Jeans study except 
that the growth of the perturbation is algebraic instead 
of being exponential. For the critical index 7 = 74/3, the 
system is stable for k > kj and unstable for k < kj. 

Some curves representing the evolution of the density 
contrast S{a) in the different cases described above are 
represented in Figs. [T][3]for illustration. 

Remark 1: for 7 > 74/3, the system is asymptotically 
stable. However, for k < fcj(O), the perturbations start 
growing before decaying. It can happen that their growth 
generates nonlinear effects that may trigger new instabil- 
ities. Therefore, their full stability is not granted. It is 
likely that the nonlinear evolution leads to a process of 
fragmentation marked by the formation of dense local- 
ized clusters in hydrostatic equilibrium (virialized struc- 
tures). These structures correspond to complete poly- 
tropic spheres that are dynamically stable since 7 > 74/3 
(see Sec. IiTbIi . 

Remark 2: it can be relevant to also treat the case 
where < (corresponding to Kj < 0). In certain 



At T = 0, it is easy to compute Sp/pb ~ <5(x, t) in sec- 
ond order perturbation theory and determine the growth 
of the skewness due to nonlinear gravitational clustering. 
This is a classical calculation in cosmology [i^ that we 
shall adapt to the present situation. For a cold over- 
damped gas, the equations of the problem are 

o e -1 

^— = SaGpb5{l + S) + —VS-\/c^, (51) 
at 



A(j) = SdGpbSa^ 



(52) 



It is convenient to measure the evolution in terms of a 
instead of t. Using Eq. (|22|) . we obtain 



a dS 
d da 



Ax = <5, X = 



SdGpbO^ 



(53) 



(54) 



In the linear approximation, the density contrast and the 
gravitational potential are given by 

(5ox,a = — —, xox,a) = — —, 55) 
a" a" 

with AC = D. To estimate the deviation to the linear 
regime, we write 



(5(x, a) = (5o(x, a)[l + e(x, a)]. 



(56) 



with (5o ^ 1 and e ^ 1. Substituting Eq. ([55]) in Eq. 
([53|. we obtain at leading order 



= + 



(57) 



Using the expressions of 5q and xo given in Eq. (j55p . the 
foregoing equation can be rewritten 



3e d fr. "^D 



(58) 
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Integrating over a, we obtain 



D 



So + ^-Vxo- (59) 



Therefore, the expression of the density contrast in sec- 
ond order perturbation theory is 



(5 = (5o + (5^ + V(5o-Vxo- 



(60) 



This equation shows that, at this order, the behavior of 
the system is non-local since the density contrast in x 
depends on the values of Sq in x' through the function 
Xo that is solution of the Poisson equation Axo = i^o- 

An interesting application of the previous result con- 
cerns the growth of the skewness of the distribution 
S{x,t). We assume that the initial distribution Si{x) is a 
random Gaussian process with zero mean and autocorre- 
lation function £,{x). Therefore (Si) = and {6f) = ^(0). 
In the linear regime, the process remains Gaussian and 
the autocorrelation increases like £,{x,a) oc In the 

nonlinear regime, the relation (5) = must be preserved 
because of mass conservation. This can be checked ex- 
plicitly at second order from Eq. Indeed, 

(S) = {So) + {S'o) + (V<5o • Vxo) = {S'o) - {S'o) = 0, (61) 

where we have used an integration by parts, and the Pois- 
son equation fM]) . to get the second equality. The inter- 
esting moment is the skewness (5^) since the initial dis- 
tribution (Gaussian) has zero skewness. From Eqs. ((56)) 
and (|59p . we obtain at lowest order 

(S^) = {Slil + 3e)) = 3{Sle) = 3(6^) + HS^VSo ■ Vxo). 

(62) 

This expression can be simplified as follows 

{S') = 3{Sl) + {\/6l ■ Vxo) = 3(^^) - {S^o) = 2('5o')- (63) 

For a Gaussian distribution, the Kurtosis is equal to 3 
i.e. {6q) — 3^(0)^. Therefore, the skewness of the density 
fluctuations is 



where, to lowest order, the variance is (5^) 
leads to our final result 



((52)2 



6. 



(64) 
C(0). This 

(65) 



B. An exact integral equation in Fourier space 

It is somewhat easier to analyze the nonlinear evolution 
of the system in Fourier space rather than in direct space. 
We shall therefore rewrite the exact equations (p4)) - ([35)) 
for the density contrast in Fourier space. We follow a 
method similar to the one developed in cosmology (49l . 

M. 



If we restrict ourselves to ideal self-gravitating Brown- 
ian particles described by an isothermal equation of state 
(jSj, the fundamental equations of the problem are 



V(5 + (1 + (5)V0 



Acf) ~ SdGpbSa^. 



(66) 



(67) 



Let us decompose the density contrast and the gravita- 
tional potential in Fourier modes 



5(x,t) = / S^{t) 



(2^ 



(x,t) 



0k(i)e 



zkx 



dk 

(2^ 



The Poisson equation (p7)) implies 



SdGpbO^ 



SdGpba'^ ' 



(68) 



(69) 



(70) 



where we recall that pf,a'' ~ dM / Sd is constant. On the 
other hand, substituting Eqs. (|68)) and ((69)) in Eq. (|66)) . 
and using Eq. (jTO)) . we obtain after simple calculations 



-k (5k + SdGpbSu 



+SdGpb / (5k'(5k-k 



kk' dk' 



k'^ (271)'' ■ 



(71) 



Symmetrizing the last term, we finally get 



SdGpb 



^k'(5k-k' 



ksT 

C(5k = 7fc^(5k + SdGpbSu 

k k' k • (k - k') 

"fc^^ Ik-k'P 



dk' 



^•(72) 



Using Eq. (|70p , we can obtain an equation for the Fourier 
components of the gravitational potential. After simple 
calculations and re-arrangement of terms, it can be writ- 
ten 



^(^k + 2(d-l)e-,^k= 

a ma 



, 2 

fc (pk 



X / '^ik+p'/'ik-p 



1 

~ 2^ 
dp 



(73) 



where we have defined p = k' — k/2. Measuring the 
evolution in terms of a instead of t, and using Eq. ()22|) . 
we obtain the equivalent equations 



/ '5k'(5k-k' 



a (i(5k 
d da 

k k' k • (k - k' 



keTk'^ 
SdGpbTna' 



-1(5] 



Ik k' 



dk' 



(74) 
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and 



da 

Hk+p-^ik-p 



2(rf-l)0k 



ksTa' 



d-2 



GMm 



k 



2GM 
dp 



(75) 



We stress that these equations are exact and closed 
while the equivalent equations in cosmology [i^ [s^ are 
not closed. These equations, however, like the original 
Snioluchowski-Poisson system ([T])-®, rely on a mean- 
field approximation. More exact (but less useful) equa- 
tions, which do not make the meanfield approximation, 
are derived in Appendix ICl for T ~ 0. 

Like in cosmology, it is possible to derive interesting 
results from these integral equations in Fourier space. We 
shall mention in particular the fc^ tail and the self-similar 
solution. 



1. Inverse cascade: the k* tail 



and (5^^^ is the next order correction. To second order, 
we get from Eq. ([71)1 with T = the equation 



a dS, 



(2) 



d da 



k • k' k • (k - k') 



|k-k'|2 



dk' 

(27r)° 



(78) 



Using Eq. ([75]) . it can be rewritten 



a dS[ 



(2) 



(2) 



2a 



2d 



i^k'i?k-k 



d da ' 

k • k' k • (k - k') 



|k-k' 



dk' 

(27r)« 



(79) 



The solution to this equation is the sum of the solution 
of the homogeneous part (which increases like 1 / a'^) and 
a particular solution (which increases like l/a^''). We 
shall only keep the particular solution which dominates 
for large times. Writing 6k = Ci^/a^'^, determining the 
constant Ck from Eq. ([79]), and reintroducing the origi- 
nal variables, we obtain the solution 



We consider the case of a cold overdamped gas (T = 
0). In the linear regime, the density contrast and the 
gravitational potential are given by 

^(-)(x,a).^, 0(-)(x,a) = ^^^g^,(76) 
where Q is solution of the Poisson equation 



AQ = D, 



Ok 



k^ ■ 



(77) 



Let us suppose that the initial power spectrum for the 
density has very little power at large scales so that 
Ps{k) (X fc" with n > 4 for small k. In the linear 
regime, the power spectrum evolves like Ps{k) cx ^"/a^'^ 
so that its shape in fc-space is not modified and is al- 
ways subdominant to fc^. However, when nonlinear ef- 
fects come into play, a "fc^ tail" develops for fc — )• (see 
below) and the contribution of the large-scale structures 
finally dominate the power spectrum. This result was 
first found in cosmology H^Hol and can be interpreted as 
a sort of "inverse cascade" like in two-dimensional turbu- 
lence. In 2D turbulence, the inverse cascade is associated 
with the formation of large-scale vortices. Similarly, in 
the present problem, it is associated with the formation 
of virialized clusters (in cosmology, these clusters corre- 
spond to galaxies and dark matter halos). It may be 
noted that the equations of 2D turbulence for the coarse- 
grained vorticity field [s^ [s^l are relatively similar to the 
Smoluchowski-Poisson system. This remark can reinforce 
the analogy between self-gravitating (Brownian) systems 
and two-dimensional turbulence |8a | . 

To see the emergence of the fc^ tail in nonlinear grav- 
itational clustering, we write 6 = (J'^-* + S'^^^ + ... where 
(5(1) = (5(^) is the density contrast in the linear theory 



r(2) 



1 



c(i)r{i) 



2 / ^k'^^Lk' 



k • k' k • (k - k') 

Ik-k'P 



dk' 



1-m 



For fc — > 0, we get 



^(2) 
''k~0 



2 



(i)|2 



dk' 

(27r)° 



a 



(81) 



This implies that the power spectrum for the density be- 
haves like Ps{k) (X k'^/a^'^ for fc -> 0. 

Of course, the same result can be obtained from Eq. 
(|75]) . Writing (f> ^ (f)^^'^ + ^f^' + ... and using Eq. ((76l) . 
we obtain 



d_g^ 
da 



+ 2(d-l)0i,') 



a 



d-2 



2GM 

.d\2 



,(L) ,(L) s rfp 

^'ik+p^U-pG(k,p)(^ 



{SdGp.a'^f f dp 
2GMa3'^-2 J ^|k+pyik-pt-(k,P) 



(82) 



Proceeding as before, the solution to this equation is 

„d-2 



2dGM 
For fc ^ 0, we get 



,(L) ,(L) ^ dp 

%l^^AL^G{k,p)^^. 



^(2) 



2dGM 



L)|2p2 



dp 

(2^ 



1 



cx 



(83) 



(84) 



which is independent on k. Using Eqs. (j70|) and 
the power spectrum for the density behaves like Ps{k) 
kia^(d~2)p^^^^ (x k^/a^"^ for fc ^ 0. 
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2. Self-similar solution in the fully nonlinear regime 



C. Spherical approximation 



At T = 0, Eq. ([74)) admits self-similar solutions of the 
form Si^{a) — S{a)Dif^. This self-similar solution is ex- 
pected to be reached for sufficiently large times when 
the initial condition becomes irrelevant. This implies 
that nonlinear gravitational clustering leads to a univer- 
sal power spectrum at late times. Again, this result is 
similar to the one obtained in cosmology |50|. However, 
in the present case, the self-similar solution is exact while 
in cosmology it is approximate (although very accurate) 
since the Zeldovich approximation must be used to close 
the integral equation for the density contrast. 

Substituting the ansatz (5k(a) — '5(a)Z?k in Eq. (|74p. 
we obtain the separate equations 



a dS 
d da 

^k'i?k-k' 



k • k' k • (k - k') 



(85) 



k' 



|k-k' 



dk' 

(27r)" 



(86) 



where fi is an arbitrary constant. Eq. ()85p can be easily 
solved yielding 

We note that S{a) — > 1/fi for large times (a — )• 0) showing 
that nonlinear effects finally stabilize the system. Equa- 
tions (|86p and (j87p provide an exact solution of the non- 
linear problem. For k — s- 0, Eq. ([55)1 reduces to 



(2^)'^' 



implying that the power spectrum for the density behaves 
like Ps{k) (X S{a)^k'^ for fc ^ 0. 

Of course, we can do the same analysis from Eq. ([75]). 
Substituting (/)k(a) = F{a)SdGpba'^Q\^ in Eq. ([75|. we 
get the separate equations 



dp 

+ 2{d - 1)F = dfia'^-^F^, 

da 



t^Qk = lJ Qik+pOik-pG(k,p)-^^, 



(89) 



(90) 



with G(k,p) = (A:/2)2 



P 



2(k-p/fc)2. This is con- 



sistent with the preceding formulae if we set S{a) ~ 
F{a)a'^-'^ and = -fc^Qk- For fc ^ 0, Eq. ^ re- 
duces to 



MQk.o = ^/ IOplV(^, 



(91) 



implying that the power spectrum for the density behaves 
like Ps{k) - k^a^^'^-^^P4,{k) (X F(a)2a2(''-2)fc4 for fc ^ 0. 

Remark: in Appendix |D1 wc show that the self-similar 
solution can be directly obtained in physical space (in- 
stead of Fourier space). 



In the nonlinear regime, it is not possible to analyt- 
ically solve Eqs. in full generality. We can, 
however, obtain a particular solution if we assume that, 
locally, the density is homogeneous. Therefore, we shall 
describe the evolution of a constant overdensity com- 
pared with the background. This is the equivalent of 
the spherical top hat in cosmology [i^, [s^l . If we assume 
(5(x,t) = 5{t), Eq. (iMl) becomes 



dt 



(92) 



Let us assume that this equation describes the evolution 
of a spherical region of mass M. and radius R{t). In that 
case, the density contrast can be written 



p M / a\d 
1 + S= — ^ 
Pb 



M \R 



(93) 



Substituting Eq. in Eq. 

we obtain 



GM 



dR _ 
^~dt ^ Rd 



or 



and using Eq. 



(94) 



This corresponds to the collapse of a uniform sphere de- 
scribed by the BSP system ([I])-©. 

This observation leads to the following scenario. If we 
start from a uniform sphere of mass M and radius oq, 
this sphere is unstable and collapses under its own grav- 
ity. The mean flow develops a finite time singularity: the 
average density Pb{t) becomes infinite in a finite time t^, 
while the radius a(t) of the sphere tends to zero. How- 
ever, this sphere is unstable to small perturbations and 
a process of fragmentation follows. The details on this 
process depend on the equation of state, e.g. on the poly- 
tropic index. If 7 > 74/3, the system is asymptotically 
stable. However, the perturbations whose wavelengths A 
are larger than the initial Jeans wavelength Aj(0) grow 
(see Sec. IIII C|) . This initial growth can trigger nonlinear 
effects and result in the formation of overdense regions. 
At first, these regions experience a cold collapse and in- 
tensify. Then, when pressure effects become important at 
high densities and small scales, localized clusters in virial 
equilibrium form. These individual clusters, which cor- 
respond to complete polytropic spheres, are stable since 
7 > 74/3- If 7 < 74/3, the collapsing sphere of radius 
a{t) is unstable to small perturbations whatever their 
wavelength. In the linear regime, the density contrast is 
small ((5^1) and the perturbations grow as described in 
Sec. IIII CI In the nonlinear regime, the density contrast 
becomes large {S ^ 1) and the system forms overdense 
regions. They correspond to localized clusters. How- 
ever, these individual clusters are unstable since 7 < 74/3 
and they experience gravitational collapse. They can col- 
lapse as uniforms spheres and fragmentate into smaller 
pieces as described in this paper. This leads to a process 
of hierarchical clustering in which clusters fragmentate 
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into smaller clusters which themselves fragmentate into 
smaller clusters and so on. Alternatively, the clusters 
can undergo a self-similar spherically symmetric collapse 
up to a Dirac peak as described in [lol - [T^ . [33 | . In that 
case, there is no fragmentation. The dynamics of these 
clusters is also interesting since they interact with each 
other. At large scales, their internal structure can be 
neglected and the problem is reduced to the dynamical 
evolution of a system of Af clusters in interaction. When 
two clusters come into contact, they merge so that their 
mass M{t) increases while their number J^{t) decreases. 
Ultimately, only one cluster remains: a Dirac peak con- 
taining all the mass (statistical equilibrium state in the 
canonical ensemble). This is reminiscent of a coarsen- 
ing process in statistical mechanics [66j . This also shares 
analogies with the process of decaying two-dimensional 
turbulence 67[. Starting from an incoherent initial con- 



mr,t)\^dr = l. 



(98) 



dition, large-scale vortices spontaneously emerge in a 2D 
incompressible flow and rapidly dominate the dynamics. 
These vortices evolve under the effect of their mutual ad- 
vection (due to long-range interactions) punctuated by 
merging events. Their number decreases as a power law 
N{t) - t-^ with an exponent ^ ~ 0.6 - 1 jgl]. It would 
be interesting to further develop the analogy between 2D 
decaying turbulence and self-gravitating Brownian parti- 
cles. 



V. QUANTUM SMOLUCHOWSKI-POISSON 
SYSTEM 

We now introduce and study a more general Smolu- 
chowski equation taking quantum effects into account. 
To that purpose, we consider a dissipative gas of 
self-gravitating Bose-Einstein condensates. Using the 
Madelung transformation, we show that the dissipative 
Gross-Pitaevskii equation is equivalent to the damped 
barotropic Euler equations with an additional quantum 
potential. In a strong friction limit, this yields the quan- 
tum barotropic Smoluchowski equation. We extend our 
previous analysis to this more general equation. 



A. Dissipative Gross-Pitaevskii-Poisson system 

We consider an interacting gas of Bose-Einstein con- 
densates at T = described by the dissipative mean field 
Gross-Pitaevskii (GP) equation [69l[70|: 

ih^ (r, t) = -—A^{r, t) + m'^tot (r, i)^(r, t) 
at 2m 

*tot(r,i) = I p{r\t)u{\r~v'\)dr', (96) 
p(r,i) = 7Vm|^(r,i)P, (97) 



Equation ([M)) is the normalization condition, Eq. ([97)1 
gives the density of the BEG, Eq. ([U^l determines the 
associated potential and Eq. ([M]) determines the wave 
function. We assume that the potential of interaction 
can be written as u = + usr where u^fi refers to 
the long-range gravitational interaction and usr to the 
short-range interaction. We assume that the short-range 
interaction corresponds to binary collisions that can be 
modeled by the effective potential usfl(r— r') = g6{r—r'), 
where the coupling constant (or pseudo-potential) g is 
related to the s-wave scattering length as through g = 
Anasfi^ /nn? (in d = 2>) [7l|. For the sake of generality, 
we allow as to be positive or negative (a^. > corre- 
sponds to short-range repulsions and Os < corresponds 
to short-range attractions). Under these conditions, the 
total potential can be written ^tot = ^ + h{p) where <& 
is the gravitational potential and 



Kp) = 9P = gNmlipl'^' 



(99) 



is an effective potential modeling short-range interac- 
tions. For the sake of generality, we will consider an 
arbitrary function h{p). We shall also assume that the 
particles are subjected to an external potential ^ext{^)- 
For illustration, we may consider a harmonic potential 
^ext = Wor'^/2 which can play the role of a trap {ujq > 0) 
or a rotation (wq = —fl^ < 0). Regrouping these re- 
sults, we obtain the dissipative Gross-Pitaevskii-Poisson 
(GPP) system 

ih-^ (r, t) = -—A^Pir, t) + m[$(r, t) 



+h{p) + $e^t(r)]^/'(r, t) + ^ In 



A$ = SdGp = SdGNm\ij\'^ 



^, (100) 



(101) 



Instead of the gravitational potential, we could consider 
an arbitrary binary potential with long-range interac- 
tions of the form Uifl(|r — r'|) in which case the Pois- 
son equation would be replaced by $ = u^jf * p where * 
denotes the product of convolution. 

We have introduced a source of dissipation in the 
Gross-Pitaevskii equation measured by the coefficient ^ 
(the ordinary Gross-Pitaevskii equation is recovered for 
^ = 0). Initially, we "guessed" the form of this dissipa- 
tion term in order to obtain, after making the Madelung 
transformation, the damped Euler equations (as we shall 
see, the parameter ^ can be interpreted as a friction coef- 
ficient in the Euler equation (jllip ). Then, we found some 
confidence in the fact that the same term was previously 
derived by Kostin [72| from the Heisenberg-Langevin 
equation describing a Brownian particle interacting with 
a thermal bath environment [93j . We can therefore con- 
sider that Eqs. p00p - (|10ip have some physical founda- 
tion. Recently, some authors |51| have proposed that 
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dark matter in the universe could be a self-gravitating 
Bose-Einstein condensate described by the GPP system 
(fT00ll - (fT0T|) with ^ = (see ^ for a short history). Here, 
we consider a more general model where damping effects 
are taken into account. In addition to cosmology, this 
model can probably find applications in other domains 
of physics. 



B. Madelung transformation 

Let use the Madelung [t^ transformation to rewrite 
the dissipative GPP system in the form of hydrodynamic 
equations. We write the wavefunction in the form 



(102) 



where A{r, t) and S{r, t) are real functions. We clearly 
have 



A 



2i VV* 



(103) 



Note that the dissipative term in the GP equation (|100p 
can be written ^S'V'. Substituting Eq. (fT02|) in Eq. (fTOOj) 
and separating real and imaginary parts, we obtain 



dS_ 
'dt 



2m 



2m~A 



C5' = 0, (104) 



(105) 



Following Madelung, we introduce the density and the 
velocity fields 



generalized Bernouilli equation for a potential flow. Tak- 
ing the gradient of Eq. (|108|) . and using the well-known 
identity (u • V)u — V(u^/2) — u x (V x u) which reduces 
to (u- V)u = V(u^/2) for an irrotational flow, we obtain 
an equation similar to the damped Euler equation with 
an additional quantum potential 

1 

^ + (u-V)u = -Vh-V^-V^e^t VQ-Cu. (110) 

ot m 

This equation shows that the effective potential h ap- 
pearing in the GP equation can be interpreted as an en- 
thalpy in the hydrodynamic equations. We can rewrite 
Eq. ([TTU)) in the form 

^ + (u-V)u = --Vp-V$-V$e.t VQ-eu, (111) 

ot p m 

where p{r,t) is a pressure. Since h{r,t) = h[p{r,t)], the 
pressure p{r,t) = p[p{r,t)] is a function of the density 
(the flow is barotropic). The equation of state p{p) is 
determined by the potential h(p) through the relation 



h'ip) 



P'ip) 



(112) 



This yields p{p) = ph{p) — H{p) where H is a primitive of 
h. In conclusion, the dissipative GPP system is equivalent 
to the hydrodynamic equations 



0, 



(113) 



^ + (u-V)u = --Vp-V$-V$e.t VQ-eu, (114) 

ot p m 



p = NmA^ ^ Nm\^P\^, u = — V5. (106) 

m 

We flrst note that the flow is irrotational since V x u = 
0. With these variables, Eqs. (fT04l) and (fTOS]) can be 
rewritten 



A$ = SdGp. 



(115) 



For a harmonic potential, W^ext = WqT. We shall refer 
to these equations as the quantum damped barotropic 
Euler-Poisson system. We note the identity 



|j+V.(pu) = 0, 



(107) 



m P 



(116) 



dS 1 

-^ + 7^{^S)^+m{^+h{p) + ^,,t) + QHS = 0, (108) 
ot Im 

where 



where P^- is the quantum stress (or pressure) tensor 



Am 



rpdtdj hip, 



(117) 



Q 



2m y/p 



Am 



Ap 
P 



1 {Vpf 



(109) 



is the quantum potential. The first equation is similar to 
the equation of continuity in hydrodynamics. The second 
equation has a form similar to the classical Hamilton- 
Jacobi equation with an additional quantum term and 
a source of dissipation. It can also be interpreted as a 



(118) 



This shows that the quantum potential is equivalent to 
an anisotropic pressure. In the classical limit ?i — )■ 0, the 
quantum potential disappears and we recover the ordi- 
nary damped barotropic Euler-Poisson system. 
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In the strong friction limit ^ — >■ +00, we can formally 
neglect the inertia of the particles in Eq. (|114[) and obtain 



The corresponding GP equation can be written 



~ --\/p - V$ - V$ext - — VQ. 
p m 



(119) 



Substituting this relation in the continuity equation 
(|113p . we finally arrive at the quantum barotropic 
Smoluchowski-Poisson system 

^1^^^-(^P + P^'^ + P^^--t + -^Q) ' (120) 
ot \ m J 



A$ = SdGp. 



(121) 



In the classical limit ?i — > 0, the quantum poten- 
tial disappears and we recover the classical barotropic 
Smoluchowski-Poisson system (IT])-©. We note that, due 
to the complex nature of the wave function, it is not pos- 
sible to take the strong friction limit directly in the dis- 
sipative GP equation (jlOOp . We must necessarily split 
this equation into its real and imaginary parts, which is 
achieved by means of the Madelung transformation, and 
take the limit ^ — -f 00 in the Euler equation (|114p cor- 
responding to the real part of the GP equation. 

Remark: a system of classical Brownian particles with 
long and short range interactions [s^ is described by the 
generalized Smoluchowski equation 



-Vp + \7pex + pV* + pV$e 



(122) 



where pid = pkBT/m is the ideal kinetic pressure and pe^ 
is the excess pressure taking into account short-range in- 
teractions. The total pressure is p = pid +Pex- In the 
present case, T = 0, so that the pressure in Eq. (|120p 
corresponds to Pex- Contrary to the kinetic pressure pid, 
the pressure pex arising from short-range interactions can 
be positive or negative (see Sec. IV C[) . In the DDFT, 
this pressure is related to the excess free energy Fex[p] 
by ypex = P^SFex/Sp d^. In the present case, accord- 
ing to Eq. (|112p . we have Vpex = p^h. Therefore, the 
excess free energy is Ff,x = J H{p) dr where 77 is a prim- 
itive of h. Similarly, the quantum pressure in Eq. (jl20p 
can be obtained from the quantum excess free energy 
pQ = (1/m) / pQdr via (p/m)VQ ^ p\76Fg/6p (see 
Appendix [E|. 



C. Equation of state 

The effective potential h{p), which takes into account 
short-range interactions (collisions) between particles, 
determines a barotropic equation of state p{p) through 
the relation (|112p . Inversely, specifying an equation of 
state p = p{p), we can obtain the corresponding function 
h{p). Let us give specific examples. 

The isothermal equation of state p ~ phgT^f / /m leads 
to an effective potential of the form 



Hp) = In p. 



m 



(123) 



2m 



AiP + {m^ + 2kBnffln\tlj\)tlj. (124) 



Interestingly, we note that a nonlinear Schrodinger equa- 
tion with a logarithmic potential similar to Eq. (|124p has 
been introduced long ago by Bialynicki & Mycielski [t^I 
as a possible generalization of the Schrodinger equation 
in quantum mechanics. 

The polytropie equation of state p = Kp'^ with 7 = 
1 -f 1/n leads to an effective potential of the form 



^7-1 



(125) 



The corresponding GP equation can be written 



^- — Ai' + m{<i> + K\^p\''^'')^, (126) 
ot 2m 

where n — K{n + l){Nm)^^'" . This is the usual gen- 
eral form of the GP equation considered in the literature 
[TSj . We may recall that classical and ultra- relativistic 
fermion stars at T = are equivalent to polytropes with 
index n = 11,^/2 = d/2 and n ~ n'^ = d respectively [s^ . 
The GPP system (|T^ - POT|) with n = 3/2 in d = 3, 
describing self-gravitating fermions beyond the Thomas- 
Fermi approximation (i.e. with the quantum potential 
retained), has been studied by Bilic et al. [t^ in relation 
with the formation of white dwarf stars by gravitational 
collapse. 

The original Gross-Pitaevskii equation corresponds to 
a potential of the form 



Kp) = 9P = 



(127) 



where is the s-scattering length (in d = 3) . This leads 
to an equation of state 



1 2 2nash^ 2 

p ^gp = — 5— P ■ 

2 m"^ 



(128) 



This is a polyropic equation of state of the form with 
n = 1, 7 = 2 and A' = g/2 = 2TrasH/m^. The GPP sys- 
tem p26p-(|1 0i p with n = 1 has been studied by Bohmer 
and Chavanis I77h79|| in the context of dark 



51 



& Harko 
matter. 

Remark: as we have previously indicated, for a Bose- 
Einstein condensate, the pressure can be negative. This 
is the case, in particular, for a BEG with quartic self- 
interaction {p Gc p^ (X IV'I^) described by the equation of 
state (|128p when the scattering length is negative. In 
terrestrial BEG experiments, some atoms like ^Li have 
a negative scattering length [Til, [s^l ■ Similarly, the con- 
stant Te/ / appearing in Eq. (|123p is just an "effective" 
temperature since it arises from a particular form of self- 
interaction and has nothing to do with the thermody- 
namical temperature (which here is T = 0). This effec- 
tive temperature can be positive or negative. 
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D. Time-independent dissipative GP equation 



If we consider a wavefunction of the form 



(129) 



we obtain the time-independent dissipative GP equation 



2m 



AV'(r) + m($(r) + h{p) + $ext(r))i/'(r) = EiP{y), 



(130) 

where -0(r) = Air) is real and p(r) = NmilP'{Y). Dividing 
Eq. ([T30| by i}j{v), we get 



m<i> + mh{p) + m<i>e2..t(r) - 



or, equivalently, 



2m 



i?, (131) 



TO(i> + mh{p) + m$ea.t + Q = 



(132) 



This relation can also be derived from the dissipative 
quantum Hamilton- Jacobi equation (jlOSp by setting S ~ 
-{E/^){1 - e"^*). Alternatively, it can be deduced from 
the quantum barotropic Euler equation (jllip which is 
equivalent to the GP equation. The steady state of the 
quantum barotropic Euler equation (jllip . obtained by 
taking dt = and u = 0, satisfies 



2™2 y ^ J 



(133) 



This generalizes the usual condition of hydrostatic equi- 
librium by incorporating the contribution of the quan- 
tum potential. Equation (|133p describes the balance be- 
tween the gravitational attraction, the repulsion due to 
the quantum potential (Heisenberg uncertainty princi- 
ple) and the repulsion (for > 0) or the attraction (for 
as < 0) due to the short-range interactions (scattering). 
This equation is equivalent to Eq. (|131|) . Indeed, inte- 
grating Eq. (fT33l) using Eq. (fm]) . we obtain Eq. (fT3T|) 
where the eigenenergy E appears as a constant of inte- 
gration. On the other hand, combining Eq. (jl33p with 
the Poisson equation (|10ip . we obtain the fundamental 
equation of hydrostatic equilibrium with quantum effects 

For a harmonic potential, A^^xt ~ cIloq. When quantum 
effects are ignored (classical limit or Thomas-Fermi ap- 
proximation), we recover the classical fundamental equa- 
tion of hydrostatic equilibrium 



(135) 



Finally, for an equation of state of the form (|128p . Eq. 
(|134p becomes 



2m^ \ y/p J 



SdGp + A$, 



(136) 



This equation has been studied in [77|, |78| in d = 3 when 
^ext = 0- There are two important limits: (i) the non- 
interacting limit corresponding to g = = (or more 
generally p = 0) and (ii) the Thomas-Fermi (TF) limit 
obtained by neglecting the quantum potential Q. 



E. Equation for the evolution of the density 
contrast 

We can now generalize the results of the first part of 
the paper to the quantum BSP system (|120p - (|12ip . In 
this paper, we assume ^ext = 0. The background flow 
(collapsing uniform sphere) studied in Sec. Ill Gl is not 
modified since the quantum potential enters (jl20p only 
through its gradient. Then, working in the comoving 
frame and extending the calculations of Sec. IIIIl we find 
that Eqs. (1M])-(|3S]) are replaced by 



.95 1 „ 
ot 

-2 



c^V(5 + (1 + (5)V(/> 



r(l + (5) r AVTTs 

2m2a2 yrri 



A(j> = SdGptSa^ 



(137) 



(138) 



where = p'{p) = p'(p6(l + S)). If we consider small 
perturbations (5 <C 1, </> ^ 1 around the homogeneous 
flow (5 = = 0; we obtain the linearized equation 



dt a2 



AS + SdGpbS 



Am?a' 



-A^5, (139) 



where = p'{pb). Expanding the solution in Fourier 
modes of the form (5(x, i) = 5\^{t)e'^^'^ ^ we obtain the 
equation 



^5 



SdGpb U = 0, (140) 



for the evolution of the density contrast in the linear 
regime. 

F. Jeans type analysis in a static frame 

Let us first make the Jeans swindle and assume that 
the above equation is valid in a static frame a = 1. In 
that case, the solution of Eq. (jl40p is (5 oc e~*"* leading 
to the dispersion relation 



\k'^ — SdGp. 



(141) 



The perturbation evolves like e^^ with an exponential 
rate 7 = —iuj. This relation, like the original quantum 
Smoluchowski equation (|120p . clearly shows the compe- 
tition between the attractive gravitational force, the (re- 
pulsive or attractive) pressure arising from short-range 



16 



interactions and the repulsive quantum pressure arising 
from the Heisenberg uncertainty principle. For the equa- 
tion of state (|128|) . the square of the velocity of sound 
is 



9P ■ 



and the dispersion relation can be rewritten 



Am? 



gpk - SdGp. 



(142) 



(143) 



In the non-interacting case g = ~ Q, the pressure 
is zero ij) = 0). The particles interact via gravity and 
they experience the effect of the quantum potential. The 
dispersion relation reduces to 



2l4 



Am? 



- SdGp, 



leading to the quantum Jeans wavenumber 



kn = 



(144) 



(145) 



This characteristic wavenumber arises due to the inter- 
play between gravity and quantum effects (Heisenberg's 
uncertainty principle). For k < kq, the system is un- 
stable and the perturbations grow with a growth rate 
7 > 0; for k > kq the system is stable and the pertur- 
bations decay with a damping rate 7 < 0. Note that the 
growth rate is maximum for k — Q (corresponding to an 
infinite wavelength A = 27r/A: — > -^00) and its value is 

Imax = SdGp/^. 

In the Thomas-Fermi approximation in which the 
quantum potential can be neglected, the particles inter- 
act via gravity and they experience a pressure due to 
short-range interactions. The dispersion relation reduces 
to 



i^uj = gpk^ — SdGp. 



(146) 



This is the classical Jeans dispersion relation of Sec. IIIBI 
For g < 0, the system is always unstable. For g > 0, the 
Jeans wavenumber is 



SdGp 



1/2 



SdG 
9 



1/2 



(147) 



This characteristic wavenumber arises due to the inter- 
play between gravity and repulsive scattering. We note 
that the Jeans wavenumber is independent on the den- 
sity (in d = 3, it can be written kj = (Gm^/ag?!^)^/^). 
For k < kj, the system is unstable and the perturbation 
grows with a growth rate 7 > 0; for k > kj the system is 
stable and the perturbation decays with a damping rate 
7 < 0. The maximum growth rate corresponds to /c = 
leading to jmax = SdGp/S,. 



In the non-gravitational case {G = 0), the dispersion 
relation reduces to 



Am^ 



gpk 



(148) 



For g > 0, the system is always stable and the perturba- 
tion decays with a damping rate 7 < 0. For g < 0, the 
critical wavenumber is 



fcn = 



~1F 



Ajr?\g\p 



1/2 



(149) 



This characteristic wavenumber arises due to the in- 
terplay between quantum pressure and attractive scat- 
tering (in d = 3, ko = (167r|as|p/m)^/^). For k < 
ko, the system is unstable and the perturbation grows 
with a growth rate 7 > 0; for k > ko the system 
is stable and the perturbation decays with a damp- 
ing rate 7 < 0. The growth rate is maximum for 
k, = {2\cl\myh'^y/^ = {2\g\pm^/h'^y/^ = ko/V2 and 
its value is 7* = c^m^/Sh, = g^ p^rn^ / ^fi^ (in d = 3, 
k, ^ (87r|a,|/9/m)i/2 and 7, = IQ-k"^ a^n^ p^ / ^m^) . 
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FIG. 4: Dimensionless dispersion relation —£V = k'*-|-2qk^ — 1 
with r = 7/70, K. = k/ko and a — g/go where 70 — SdGp, 
ko = {ASdGpm'^/h'^y/* and go ^ {SdGh ^ / pm^f^'^ . The 
Jeans wavenumber is Kj = — q -I- y/oP- 
maximum growth rate is ^F, — \ + o? 



-I- 1. For g < Q, the 
reached for = —a. 
The TF limit corresponds to <ti the non-gravitational 
limit to 3> l/ja] and the non-interacting limit to |a| <g; 1. 



We now consider the general case. With the previous 
relations, the dispersion relation (|141[) can be rewritten 



k^ P_ 



(150) 



SdGp 

The pulsation vanishes at the critical Jeans wavenumber 

(151) 

with -|- when > and — when < 0. The previ- 
ous results can be recovered as particular cases of this 
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general relation. For k < kc, the system is unstable 
and the perturbation grows with a growth rate 7 > 0; 
for k > kc the system is stable and the perturbation 
decays with a damping rate 7 < 0. For > 0, 7 de- 
creases monotonically with k. Accordingly, the growth 
rate is maximum for fc = (infinite wavelength) leading 
to jmax ~ SdGp/S,. For < 0, 7 achieves a maximum 
value 



SdGp 



at 



4fc4 



1/2 



SdGp 



1 4 



2|fc}l 



2\cl\m' 



1/2 



(152) 



(153) 



The exponential rate 7 = — zw is plotted as a function 
of the wavenumber k in Fig. 2] for repulsive {g > 0) and 
attractive {g < 0) self-interactions. 

Remark: if we perform the Jeans stability analysis on 
the quantum damped barotropic Euler-Poisson system 
(|113p - pi5p . we find that the dispersion relation is given 
by Eq. (|141[) where i£,uj is replaced by uj{uj + The 
discussion is then similar to the one given in Ref. [ssf . 



G. Jeans type analysis in a collapsing frame 



We now return to Eq. (|140p and study the growth 
of perturbations in the collapsing frame. Taking Cs — 
in Eq. (|140p . we can define a quantum time-dependent 
Jeans wavenumber 



I^ASdGpbrn^ 



1/4 



(154) 



Recalling Eq. ([TT]). we can write kg = /cqo^^^''^/^ with 
KQ = {iSdGpba'^m'^ /h^y/"^ . The quantum Jeans length 
Xq = 27r/fcj behaves like a"'^''" For d > 4 (resp. d < 
4), it decreases (resp. increases) with time so that the 
system becomes unstable at smaller and smaller (larger 
and larger) scales as the cloud collapses. For d = 4, the 
quantum Jeans length is constant: kj = nj. 

Measuring the evolution in terms of a instead of t and 
using Eq. (f22|) . we can rewrite Eq. (|140p in the form 



2 1,4 



dS 



da a \4:SdGpbm^a'^ SdGpta 



- 1 U = 0. (155) 



For short times, writing a = ao(l ^ c) with e ^ 1, we 
obtain 



'M^i-d( 



S{ao) 



\kQ{o)^ kj{oy 



1 e, (156) 



where fcg(O) and fcj(O) are the quantum and classical 
Jeans wavcnumbcrs (|154p and (|^T|) at t = 0. For short 



times, the perturbations start to grow if fc < fcc(O) and 
start to decay if fc > fcc(O) where kc is defined by Eq. 
(|15ip . This result is valid for any equation of state. To 
describe larger times, it is necessary to specify the equa- 
tion of state. Assuming that p{p) is given by a polytropic 
equation of state and introducing the notations defined 
previously, Eq. (jl55p becomes 




Let us first consider the non- interacting case Cs ~ 0- 
Eq. (|157p reduces to 



dS d fc4 



1 U = 0. 



If (i 7^ 4, the solution of this equation is 



-I dk 

o(a) (X —re « 



(158) 



(159) 



For short times, we recover Eq. (jl56p . This equation 
clearly shows the initial effect of the Jeans wavelength. 
Indeed, the perturbation starts to grow if fc < kgiO) and 
starts to decay if fc > fcc3(0). This result is independent 
on the value of the dimension d. On the other hand, for 
large times, the effect of the quantum Jeans scale disap- 
pears and the evolution of the system is only controlled 
by the value of d. For c? < 4, we find that (5(a) — >■ 
very rapidly as a 0. The system becomes asymptoti- 
cally stable to all wavelengths. For 0? > 4, we find that 
d{a) oc a~'^ — )■ -I-cxd as a — )■ 0. The system becomes 
asymptotically unstable to all wavelengths and behaves 
like in a cold classical gas. We note that quantum me- 
chanics has a stabilizing role with respect to gravitational 
collapse when d < A. However, it is not able to prevent 
gravitational collapse when c? > 4. Therefore, the di- 
mension d = 4 is critical. We have already reached this 
conclusion in [8l|, HH for fermion stars. For d = 4, the 
solution of Eq. (|158p is 



S{a) cx a^'^*/^ 



(160) 



For fc < Kg, the perturbation increases in time and the 
system is unstable (for fc <ti kq, the perturbation grows 
like in a cold classical gas). For k > kq, the perturbation 
decreases in time and the system is stable. 

Comparing Eqs. p58p and we note that the de- 
pendence in a is the same in the two equations when 
7 75/3 = (d + 2)/(i, i.e. n = = d/2. This is the 
index of a classical fermionic gas |8ll.l82j . Concerning the 
dependence in a, a cold bosonic gas behaves similarly to 
a cold fermionic gas. However, the dependence in fc is 
different (fc^ instead of fc) as well as the expression of 
the Jeans length (fcg instead of fcj). We also note that 
75/3 ~ 74/3 = 3/2 for d = 4, which is the critical di- 
mension of a self-gravitating quantum gas [8l|, [s^l ■ For 
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7 = 75/3 1 Eq. (|157p reduces to 

If d ^ 4, the solution of this equation is 

(5(a) oc -re (^-'')»"-'' . (162) 
a" 

For short times, we recover Eq. (|156p . The perturba- 
tion starts to grow if < fcc(O) and starts to decay if 
k > kc{0). For large times, the system is asymptoticaUy 
stable for d < 4 and asymptotically unstable for d > A. 
For d = 4, the solution of Eq. (|16ip is 

S{a) oc a (163) 

The system is stable for k > k^ and unstable for k < k^- 
Let us finally treat the general case. We first assume 
d 7^ 4. For 7 7^ 74/3, we get 

S(a) oc -;e -h^'-^^-'-\''i(^(^-^)-^-') , (164) 
and for 7 = 74/3, we obtain 

5(a) oc -ie"'=«<^"'"*-'a<''^'/'''-i). (165) 
a" 

We now assume d = 4. For 7 7^ 74/3 = 3/2, we get 

4fc2^6-47 

(5(a) (X a<'='/'"Q-^)e^^^<^^, (166) 

and for 7 = 74/3 = 3/2, we obtain Eq. (|163p . 

The short time evolution (a — > ao) is given by Eq. 
(|156p . The perturbation starts to grow if fc < fcc(O) 
and starts to decay if fc > fcc(O). To analyze the evo- 
lution for large time (a — > 0), it is convenient to start 
directly from the differential equation (|157p . (i) Let us 
first assume d < A. For 7 > 75/3, (5(a) like in Eq. 
(|47)) and for 7 < 75/3, (5(a) ^ like in Eq. (|T59)) . For 
7 = 75/3, 5(a) — > like in Eq. (|162p . In all cases, we 
find that 5 — >■ very rapidly so that the system is al- 
ways asymptotically stable. The stabilization is due to 
the polytropic pressure for 7 > 75/3, to the quantum 
pressure for 7 < 75/3 and to both for 7 = 75/3. (ii) Let 
us consider d > 4. In that case, the quantum effects 
become asymptotically negligible and we are led back to 
the study of Sec. IIII CI For 7 7^ 74/3, the asymptotic 
evolution of the density contrast is given by Eq. (|T7l). 
For 7 > 74/3, (5(a) — >■ very rapidly and the system is 
stable; for 7 < 74/3, (5(a) oc a^'^ — > +00 and the system 
is unstable. For 7 = 74/3, the evolution of the density 
contrast is given by Eq. (|50p . The system is stable for 
k > kj and unstable for k < kj. (iii) Let us finally as- 
sume d = A. For 7 > 3/2, we find that (5(a) — !• like in 
Eq. (|47l) so the system is stable. The quantum pressure 



is asymptotically negligible so the stabilization is due to 
the polytropic pressure. For 7 < 3/2, the polytropic 
pressure is asymptotically negligible. We find that (5(a) 
behaves like in Eq. (|160p so that the system is stable for 
k > kg and unstable for k < kg. The stabilization is 
due to the quantum pressure. Finally, for 7 = 3/2, the 
density contrast behaves like in Eq. (|163p so that the 
system is stable for k > kc and unstable for k < kc- The 
stabilization is due to the quantum pressure and to the 
polytropic pressure. 

Remark 1: a standard self-gravitating BEC with 7 = 2 
is always stable since 2 > 74/3 in any dimension of space. 
In that case, kj = {SdG / gY^'^ . 

Remark 2: the previous discussion assumes that > 
0. If the self-interaction is attractive (c^ < 0, Kj < 0), 
the system becomes unstable in the following cases: (i) 
if (i < 4 and 7 > 75/3; (ii) if (i < 4 and 7 = 75/3 and 
k < (^y 141)1/2; (j-J if d> 4; (iv) if d = 4and7 > 3/2. 
The other results are unchanged. 

VI. CONCLUSION 

In this paper, we have studied the dynamical evolu- 
tion of a collapsing cloud of classical and quantum self- 
gravitating Brownian particles. A spatially homogeneous 
sphere is unstable and collapses under its own gravity. 
This leads to a finite-time singularity in which the system 
has a vanishing radius and an infinite density. However, 
the collapsing sphere (background fiow) may be unstable 
to small perturbations and a process of fragmentation 
takes place. This leads to the formation of dense local- 
ized clusters. We have derived an exact set of equations 
([Ml) - ((55)) describing the evolution of the density contrast 
in the comoving frame. In the linear regime, the evolu- 
tion of the perturbations can be described analytically 
for any equation of state. Specific attention has been 
given to the polytropic equation of state for the sake of 
illustration. If 7 > 74/3, the system is asymptotically 
stable. However, the perturbations initially grow if their 
wavelength A is larger than the Jeans scale Aj(0). This 
initial growth can trigger nonlinear effects and result first 
in a cold collapse then in the formation of dense local- 
ized clusters in virial equilibrium when pressure effects 
become important. They correspond to stable steady 
states of the BSP system H])-© in which gravitational 
forces are balanced by pressure forces. If 7 < 74/3, the 
system is unstable. Small perturbations grow in the lin- 
ear regime whatever their wavelength. In the nonlinear 
regime, these perturbations amplify and form dense lo- 
calized clusters. Their structure and evolution can be 
described by the BSP system (P)-©, neglecting locally 
the contraction of the background flow. These clusters 
are unstable and they experience gravitational collapse. 
They can break into smaller fragments which themselves 
break into smaller fragments in a hierarchical manner. 
Alternatively, they can undergo a spherically symmet- 
ric self-similar collapse up to a Dirac peak. At large 



19 



scales, we can ignore their internal structure and we are 
led to a reduced dynamical system of Af clusters in in- 
teraction. The clusters merge when they come into con- 
tact and collapse on each other. Therefore, their mass 
increases while their number decreases. This is reminis- 
cent of a coarsening process in statistical mechanics or in 
2D decaying turbulence. We have given preliminary an- 
alytical results valid in the nonlinear regime. A more de- 
tailed description of the nonlinear regime requires numer- 
ical simulations. This is a problem left for future works. 
Throughout our study, we have found interesting analo- 
gies with cosmology (in a universe experiencing a "big- 
crunch" ) [4^, [5^ and with stellar formation (in a collaps- 
ing molecular cloud) [H,!!^!!^]. We have shown that the 
methods initially devised in astrophysics (Euler-Poisson) 
could be extended to a different context (Smoluchowski- 
Poisson) leading to interesting new results and applica- 
tions. We have also stressed the specificities of our ap- 
proach and the differences with the astrophysical stud- 
ies. Since the equations that we have derived 
for self-gravitating Brownian particles are simpler than 
those used in cosmology and astrophysics (because iner- 
tial effects are neglected in our model), these equations 
could be used as a valuable prototype to test the theories 
that have been developed in other domains to study the 
process of self-organization in complex media. We hope 
to develop these issues in future works. 



Appendix A: Repulsive interactions 

It can be of interest to consider the case of Brownian 
particles with repulsive interactions (see Ref. [Ill for fur- 
ther discussions and applications). The basic equations 
of the paper remain valid provided that G is replaced 
by — G. Considering a spatially homogeneous solution of 
the BSP system, as in Sec. Ill C[ we now find that the 
sphere expands, instead of contracting, as described in 
Ref. [3^. Then, from the results of Sec. IIII Bl properlv 
adapted to this new situation (i.e. making the substitu- 
tion G —G), we find that the homogeneous solution 
is always stable since S < according to Eq. (HDl). This 
is consistent with the observation made in Ref. [38| that 
the homogeneous expanding sphere is an "attractor" of 
the T = dynamics for large times. Indeed, it is found in 
Appendix E of [3S] that the inhomogeneous density pro- 
files p(r, t) converge towards the homogeneous expanding 
sphere for t — > +00. 



Appendix B: Cosmological constant 

If we introduce a "cosmological constant" A in the 
Poisson equation (this terminology is abusive here since 
we are not dealing with cosmology) and write A$ = 
SdGp — A, the equation ([^^ for the scale factor a{t) 



becomes 



da _ GM 



1 , 

-Aa. 

a 



(Bl) 



A static uniform sphere of density pf, — K/SdG and ra- 
dius a = {dG M / A)^ / is a particular solution of the equa- 
tions. To investigate its stability, we set a(t) = a + da{t) 
where Sa(t) is a small perturbation. Substituting this 
decomposition in Eq. (jBl[) and keeping only terms that 
are linear in da{t), we get 



dSa 
~dt 



A5a, 



(B2) 



so that the perturbation grows as 5a oc e^*/^ . Therefore, 
the static solution is strongly unstable similarly to the 
Einstein static universe in cosmology. 



Appendix C: Fourier transform of the equations of 
motion 



In this Appendix, we derive an exact equation for the 
Fourier transform of the discrete density contrast 5^^^ (t) 
directly from the equations of motion of the overdamped 
Brownian particles at T = 0. This derivation closely 
follows the derivation given by Peebles [4^ and Padman- 
abhan [soj in cosmology for particles with inertia. At 
T = 0, the equations of motion are 



^ dt 



A$rf = AnGpd, 



(CI) 



where Pd{'^,t) = "m^^Soly^ — Yi{t)\ is the discrete den- 
sity field expressed as a sum of Dirac distributions 
and $d(r,<) is the corresponding gravitational potential. 
Making the changes of variables and (P^ . we find 
that the equations of motion in the comoving frame are 



6 



dt 



j'^i<f>d, A4id = SdGpba Sd- 



(C2) 



The discrete density field can be written 



1 + Sd{x,t) = 



_ Pd ix,t) V yr^ 



y"(SD[x-XT(t,q)]dq, (C3) 



where XT(t,q) is the (Lagrangian) position at time t of 
the particle initially located at x = q (the subscript T 
stands for "trajectory"). Using Eq. jCSj, the Fourier 
transform of the density contrast is given by 



dq ^ 



(C4) 



Taking the time derivative of this expression and using 
Eqs. (jC2p and (fZO)) . we find that the equation satisfied 
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id) 



SdGpbS, 



kk' 



SdGpb 
2 

t- (k- 



/ 

k') 



Ad) Ad) 

"k' "k-k' 



Ik-k 



f\2 



rfk' 

(ITa with T 



(C5) 



This equation has the same form as Eq. (17^ with T = 0. 
However, its interpretation is different because, in the 
present derivation, we have not made the mean field ap- 
proximation. Therefore, Eq. (|C5|) is exact and contains 
the same information as Eq. (|C2p . It is vahd for the 
Fourier transform of the discrete density contrast (5d(x, t) 
expressed in terms of Dirac distributions while Eq. ((72)) 
is valid for the Fourier transform of the smooth (locally 
averaged) density contrast 5{x,t). Although exact, Eq. 
(|C5[) is not very useful in practice since it is difficult to 
deal with Dirac distributions. Indeed, it is easier to di- 
rectly integrate the equations of motion (|C2p that are 
equivalent to Eq. (|C5|) . By contrast, Eq. ([7^ relies on 
a mean field approximation but it applies to a smooth 
density field which is more physically relevant. 

The discrete density pd{r,t) and the discrete density 
contrast t) satisfy equations similar to Eqs. ([T]) and 
([M]) with T = constructed with the discrete gravita- 
tional potentials 4>d(r,t) and ^^(x, t). Although these 
equations are formally similar, their interpretation is dif- 
ferent as explained previously. This is the same distinc- 
tion as between the Klimontovich (exact) and the Vlasov 
(mean field) equations in plasma physics and stellar dy- 
namics. More details can be found in Ref. [s^. 

Finally, it is possible to include stochastic forces (when 
T 7^ 0) in the equations of motion. In the inertial frame, 
Eq. (|C1|) is replaced by 



'2fcBT 



R,(i), 



(C6) 



where Ri(t) is a white noise satisfying (Ri(i)) = and 



(Rf{t)R'^{t')) = SijSaisS^t — t'). In the comoving frame, 
Eq. (|C2[) is replaced by 




(C7) 



Following Ref. [86[, it is possible to derive exact equa- 
tions satisfied by the discrete density field and by the 
discrete density contrast. When a mean field approxi- 
mation is implemented, we find that the average density 
field and the average density contrast satisfy the equa- 
tions ([I])-© and derived in the main text (in 
the isothermal case). 



= SdGpbSa^ 



(D2) 



Measuring the evolution in terms of a rather than in 
terms of t, and using Eq. (P^. we can rewrite Eq. (jPip 
as 

~ + S + S' = -^^y5-Vc^. (D3) 

We consider a solution of the form 

<5(x, a) = (5(a)-D(x), 0(x, a) = F{a)SdGpba'^Q{^), 



with 



(D4) 



F{a)^a'-''5{a), AQ = D(x). (D5) 



Substituting this ansatz in Eq. (jDSj, we obtain the sep- 
arate equations 



-- + %)=M^^(a), 



AiD(x) = -D^(x) - VD • VQ. 



(D6) 



(D7) 



These equations are equivalent to Eqs. (1851) and 
The equation for 5{a) can be solved as in Sec. IIVB 21 
On the other hand, combining Eq. (jD7|) with the Poisson 
equation (|D5p -b. we find that I?(x) satisfies the nonlinear 
differential equation 



2V • 



AlnZ) 



D = 0. 



(D8) 



Appendix E: //-theorem 

The free energy associated with the dissipative GPP 
system (|100p - (jlOip . or equivalently with the dissipative 
quantum barotropic Euler- Poisson system (|113p - ()115p . 
can be written 

F = Q, + Qq + U + W + W,,t. (El) 

The first two terms correspond to the total kinetic energy 



e 



2m 



iVV-pdr. 



(E2) 



Using the Madelung transformation, it can be decom- 
posed into the "classical" kinetic energy 



Appendix D: Self-similar solution in physical space 

At T = 0, the equation for the density contrast is 

i— = SdGpb5{l + 6) + —V5-Vcp, (Dl) 
at 



and the "quantum" kinetic energy 



On ^ — I pQ dr. 

m 



(E3) 



(E4) 
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Substituting Eq. (|109|) in Eq. (jE4p , we obtain the equiv- 
alent expressions 



Or 



2m2 



f{Vpf 



2m? J ' ^' ' 8m2 j p 

The third term is the internal energy 



dr. (E5) 



U 



p I dpi dr 

Pi 



= J [ph{p)-p{p)] dr = J H{p)dv, (E6) 

where we have used Eq. (|112p and integrated by parts to 
obtain the second equality. For a polytropie equation of 
state p = Kp'^ , it takes the form 



U = / p^ dr ~ / pdr. 

1 -'^ J 1 J 



(E7) 



In particular, for the potential (jl27p . using Eq. ()128p . we 
get 



U =^g I p^dr= ^^°3^ / p2 dr. (E8) 



The fourth term is the gravitational potential energy 



W j p'^dr. 



(E9) 



The fifth term is the external potential energy 

Z"*-*^ (ElO) 

For a harmonic potential, Wext = (l/2)i^o-^ where / 
is the moment of inertia (jFip . It is easy to establish 
[73 that JGc = ^[v?/2)5pdr + / pu • dudr, 5Qq = 
(1/m) / QSpdr, SU = J h{p)Spdr, dW = J <t>Spdr and 
SWext = J ^extSp dr. Taking the time derivative of F, 
using the previous relations and the hydrodynamic equa- 
tions (jll3p - (|115p . we easily obtain the i/-theorem 



F = pu^ dr < 0. 



(Ell) 



This result remains valid if u is not a potential flow since 
u- (ux (V X u)) = 0. For a steady state, F = 0. Equation 
(jElip implies u = and we recover the condition of 
hydrostatic equilibrium (|133|) . 

The free energy associated with the quantum BSP sys- 
tem p20| - (fT2T1) is 



F = eQ + u + w + We: 



(E12) 



Taking the time derivative of F and using the previous 
relations, we easily obtain the i7-theorem 



F = - I ^iVp + + pV^ext + -VQf dr < 0. 
' t_p m 



(E13) 



It can also be obtained from Eq. (jEllI) , using Eq. (jll9l) . 
For a steady state, F = 0. Equation (|E13p implies that 
the term in parenthesis vanishes leading to the condition 
of hydrostatic equilibrium (jl33p . 

These results indicate that the system will converge 
towards a steady state of the dynamical equations that 
is a (local) minimum of free energy at fixed mass (the 
mass M = J pdr is a, conserved quantity). Maxima or 
saddle points of free energy are unstable. If several local 
minima of free energy exist, the selection will depend on 
a notion of basin of attraction. We are therefore led to 
considering the minimization problem 



min{F[p, u] 



M}. 



(E14) 



An extremum of free energy at fixed mass is given by the 
variational principle 6F—aSM = where a is a Lagrange 
multiplier taking into account the mass constraint. This 
gives u = and the condition 



m(f> + m^ext + mh(p) = ma. 

2m y/p 



(E15) 



This equation is equivalent to the steady state equation 
(|13ip provided that we make the identification 



a ~ E /m. 



(E16) 



This shows that the Lagrange multiplier (chemical poten- 
tial) in the constrained minimization problem associated 
with (|E14[) is equal to the eigenenergy E by unit of mass. 
On the other hand, considering the second order varia- 
tions of free energy, we find that the steady state is stable 
iff 

<52F=i j h'{p){5pfdr+]^ j SpS^dr 



8m'^ 



Sp_ 



dr 



8m2 



.3/2 



(Spf dr > 0, 



(E17) 



for all perturbations that conserve mass: J 6pdr = 0. 

Remark: for ^ = 0, the functional F becomes the total 
energy Etot- According to Eq. (jElip . the total energy 
Etot is conserved by the GPP (or quantum barotropic 
Euler-Poisson) system. In that context, the minimization 
problem (|E14[) provides a condition of nonlinear dynam- 
ical stability for a steady state of the GPP (or quantum 
barotropic Euler-Poisson) system f7l\ . 



Appendix F: Virial theorem 

In this Appendix, we establish the time-dependent 
virial theorem associated with the dissipative quantum 
barotropic Euler-Poisson system (|113p - (|115p . The mo- 
ment of inertia is 



pr' 



dr. 



(Fl) 
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Taking its time derivative and using the continuity equa- 
tion (|113p . we obtain 



1 = 2 I pr ■ udr. 



(F2) 



With the aid of the continuity equation, the quantum 
barotropic Euler equation (|114[) can be rewritten 

d 

— (pu) + V(pu®u) 
= -Vp - pV$ - pV$e.t - ^VQ - ^pu. (F3) 



Taking the time derivative of Eq. (jF2p , substituting Eq. 
(|F3[) . and using integrations by parts, we obtain the time- 
dependent virial theorem 

i/+ = 2(6^ + QQ) + d j pdv + Wu + Ve,f (F4) 

For a polytropic equation of state p = Kp'^ ^ using Eq. 
(|Ff|) . we have j pdr = (7 - l)U . To obtain Eq. (|F4| . 
we have used the following identities. First, it can be 
established after a few manipulations (using essentially 
integrations by part [t^]) that 



— r- VQdr = 2Qc 

m 



(F5) 



On the other hand, we have introduced the virial of the 
gravitational force 



pr ■ V<i> dv. 



(F6) 



It can be shown that W"\ = ('^ " 2)W^ if d 7^ 2 and 
W^i = -GNP/2 if d = 2 [3. Finally, we have introduced 
the virial of the external force 



pr ■ V$ej;t dr. 



(F7) 



For a harmonic potential, Vext = ^'-^qI = ^2VFext- 

In the strong friction limit ^ — > +00 in which u = 
0(1/^), the time-dependent virial theorem takes the form 

= 2QQ + d j pdr + W,^ + Vext. (F8) 

It can be directly obtained from the quantum BSP system 

For a steady state, / = / = and u = 0, we obtain 
the equilibrium virial theorem 



2Qc 



d / pdr + Wu + Vext 



(F9) 



On the other hand, the free energy reduces to F = 6q + 
U + W + Wext- Finally, multiplying the steady state 
equation p32p by p/m and integrating over the whole 
configuration, we obtain the general identity 



2W + J ph{p) dr + Wext + Oq = NE. 



(FIO) 



For a polytropic equation of state, using Eqs. (|125p and 
(|Ff|) . we find that / ph{p) dr = -fU . 



[1] C. Sire and P.H. Cliavanis, Collapse and evaporation of a 
canonical self-gravitating gas in Proceedings of the 12th 
Marcel Grossmann Meeting (World Scientific, Singapore, 
2010) arXiv:1003.1118 

[2] P.H. Chavanis, A&A 356, 1089 (2000) 

[3] E.F. Keller, L.A. Segel, J. Theor. Biol. 30, 225 (1971) 

[4] A. Dominguez, M. Oettel, S. Dietrich, Phys. Rev. E 82, 
011402 (2010) 

[5] G. Giacomin, J. Lebowitz, Phys. Rev. Lett. 76, 1094 
(1996) 

[6] D. Holm, V. Putkaradze, Phys. Rev. Lett. 95, 226106 
(2005) 

[7] J.H. Jeans, Phil. Trans. A 199, 49 (1902) 
[8] J.H. Jeans, Astronomy and Cosmogony (Cambridge Uni- 
versity Press, Cambridge, 1929) 
[9] J. Binney, S. Tremaine, Galactic Dynamics (Princeton 
Series in Astrophysics, 1987) 
[10] P.H. Chavanis, C. Rosier, C. Sire, Phys. Rev. E 66, 
036105 (2002) 

[11] C. Sire, P.H. Chavanis, Phys. Rev. E 66, 046133 (2002) 
[12] P.H. Chavanis, C. Sire, Phys. Rev. E 69, 016116 (2004) 
[13] C. Sire, P.H. Chavanis, Phys. Rev. E 78, 061111 (2008) 
[14] P.H. Chavanis, C. Sire, Phys. Rev. E 73, 066103 (2006) 
[15] A.S. Eddington, Internal Constitution of the Stars (Cam- 



bridge University Press, Cambridge, England, 1926) 
[16] S. Chandrasekhar, An Introduction to the Theory of Stel- 
lar Structure (Dover, New York, 1939) 
[17] R. Ebert, Z. Astrophys. 37, 217 (1955) 
[18] W.B. Bonnor, Mon. Not. R. astr. Soc. 116, 351 (1956) 
[19] W.H. McCrea, Mon. Not. R. astr. Soc. 117, 562 (1957) 
[20] W.B. Bonnor, Mon. Not. R. astr. Soc. 118, 523 (1958) 
[21] P. Ledoux, T. Walraven, Handbuch der Physik (Springer- 

Verlag, Berlin, 1958), Vol. 51, p. 353 
[22] V.A. Antonov, Vest. Leningr. Gos. Univ. 7, 135 (1962) 
[23] D. Lynden-Bell, R. Wood, Mon. Not. R. astr. Soc. 138, 
495 (1968) 

[24] S. Yabushita, Mon. Not. R. astr. Soc. 140, 109 (1968) 
[25] G. Horedt, Mon. Not. R. astr. Soc. 151, 81 (1970) 
[26] J. Katz, Mon. Not. R. astr. Soc. 183, 765 (1978) 
[27] T. Padmanabhan, Phys. Rep. 188, 285 (1990) 
[28] P.H. Chavanis, Astron. Astrophys. 381, 340 (2002) 
[29] P.H. Chavanis, Astron. Astrophys. 386, 732 (2002) 
[30] P.H. Chavanis, Astron. Astrophys. 401, 15 (2003) 
[31] S. Childress, J.K. Percus, Math. Biosci. 56, 217 (1981) 
[32] M.A. Herrero, J.J.L. Velazquez, J. Math. Biol. 35, 177 
(1996) 

[33] M.P. Brenner, P. Constantin, L.P. Kadanoflt, A. Schenkel, 
S.C. Venkataramani, Nonlinearity 12, 1071 (1999) 



23 



[34] 
[35] 
[36] 

[37] 
[38] 
[39] 
[40] 
[41] 
[42] 

[43] 

[44] 
[45] 

[46] 
[47] 

[48' 
[49' 



J C. Sire, P.H. Chavanis, Phys. Rev. E 69, 066109 (2004) 
i] L. Acedo, Europhys. Lett. 73, 5 (2006) 
;] N. Kavallaris, P. Souplet, SIAM J. Math. Anal. 40, 1852 
(2009) 

] P. Lushnikov, Physics Letters A 374, 1678 (2010) 
;] P.H. Chavanis, C. Sire, Phys. Rev. E 83, 031131 (2011) 
] C. Hunter, Astrophys. J. 136, 594 (1962) 
i] L. Mestel, Q. Jl R. astr. Soc. 6, 161 (1965) 
] M.V. Pension, Mon. Not. R. astr. Soc. 144, 425 (1969) 
l] R.B. Larson, PhD. Thesis, California Institute of Tech- 
nology (1968) 

1] J. A. Fillmore, P. Goldreich, Astrophys. J. 281, 1 (1984) 
] E. Bertschinger, Astrophys. J. 58, 39 (1985) 
] D. Lynden-BeU, J.P.S. Lemos, Mon. Not. R. astr. Soc. 

233, 197 (1988) 
] F. Hoyle, Astrophys. J. 118, 513 (1953) 
■] C. Hunter, Astrophys. J. 139, 570 (1964) 
1] W.B. Bonnor, Mon. Not. R. Astron. Soc. 117, 104 (1957) 
i] P. J.E. Peebles, The Large-Scale Structure of the Universe 
(Princeton University Press, Princeton, 1980) 
[50] T. Padmanabhan, Statistical Mechanics of Gravitating 
Systems: an Overview, in Long-Range Interacting Sys- 
tems, edited by T. Dauxois, S. Ruffo and L. Cugliandolo, 
Les Houches Summer School 2008, (Oxford: Oxford Uni- 
versity Press, 2009) 
[51] C.G. Bohmer, T. Harko, J. Cosmol. Astropart. Phys. 06 
(2007) 025 

P.H. Chavanis, Physica A 361, 81 (2006) 
P.H. Chavanis, Phys. Rev. E 68, 036108 (2003) 
P.H. Chavanis, Eur. Phys. J. B 62, 179 (2008) 
P.H. Chavanis, Physica A 390, 1546 (2011) 
P.H. Chavanis, L. Delfini, Phys. Rev. E 81, 051103 (2010) 
P.H. Chavanis, Eur. Phys. J. B 52, 433 (2006) 
P.H. Chavanis, C. Sire, Physica A 387, 4033 (2008) 
A. Einstein, Sitzungsber. Preuss. Akad. Wiss. 1, 142 
(1917) 

E.A. Spiegel, in Gravitational Screening, edited by 
A. Harvey (Springer- Verlag, Heidelberg, 1998); e-print 



[52' 
[53' 
[54 
[55' 
[56' 
[57 
[58' 
[59; 

[60; 



|arXiv:astro-ph/9801014 



[61] M. Kie sshng, Adv. Appl. Math 31, 132 (2003); see also 

e-print |arXiv:astro-ph/9910247] 
[62] A.S. Eddington, Mon. Not. R. Astron. Soc. 90, 668 

(1930) 

[63] E.R. Harrison, Rev. Mod. Phys. 39, 862 (1967) 
[64] W.H. McCrea, E.A. Milne, Quarterly J. Math. 5, 73 
(1934) 

[65] S. Weinberg, Gravitation and Cosmology (John Wiley & 

Sons, 1972) 
[66] A. Bray, Adv. Phys. 43, 357 (1995) 

[67] G.F. Carnevale, J.C. McWilhams, Y. Pomeau, J.B. 

Weiss, W.R. Young, Phys. Rev. Lett. 66, 2735 (1991) 
[68] C. Sire, P.H. Chavanis, Phys. Rev. E 61, 6644 (2000) 
[69] E.P. Gross, Ann. of Phys. 4, 57 (1958); Nuovo Cimento 

20, 454 (1961); J. Math. Phys. 4, 195 (1963) 
[70] L.P. Pitaevskii, Sov. Phys. JETP 9, 830 (1959); ibid 13, 

451 (1961) 

[71] F. Dalfovo, S. Giorgini, L.P. Pitaevskii, S. Stringari, Rev. 

Mod. Phys. 71, 463 (1999) 
[72] M.D. Kostin, J. Chem. Phys. 57, 3589 (1972) 
[73] E. Madelung, Zeit. F. Phys. 40, 322 (1927) 
[74] I. Bialynicki-Birula, J. Mycielski, Ann. Phys. 100, 62 

(1976) 

[75] C. Sulem, P.L. Sulem, The Nonlinear Schrodinger Equa- 



[76] 

[77] 
[78] 
[79] 
[80] 

[81] 
[82] 
[83] 
[84] 

[85] 



[87] 



[89] 
[90] 



[9i: 



[92] 



[93] 



tion (Springer, 1999) 

N. Bilic, R.J. Lindebaum, G.B. Tupper, R.D. Viollier, 
Phys. Lett. B 5 15, 105 (2001) 

P.H. Chavanis, [arXiv:1103.2050j 

P.H. Chavanis, L. Delfini, |arXiv:1103.2054] 
P.H. Chavanis, arXiv:1103.2698 

P.O. Fedichev, Yu. Kagan, G.V. Shlyapnikov, J.T.M. 

Walraven, Phys. Rev. Lett. 77, 2913 (1996) 

P.H. Chavanis, Phys. Rev. E 69, 066126 (2004) 

P.H. Chavanis, Phys. Rev. D 76, 023004 (2007) 

R. Robert, J. Sommeria, Phys. Rev. Lett. 69, 2776 (1992) 

P.H. Chavanis, J. Sommeria, R. Robert, Astrophys. J. 

471, 385 (1996) 

P.H. Chavanis, Statistical Mechanics of Two-Dimensional 
Vortices and Stellar Systems, in Dynamics and Ther- 
modynamics of Systems with Long-Range Interactions, 
edited by T. Dauxois, S. Ruffo, E. Arimondo and 
M. Wilkens, Lectures Notes in Physics 602 (Berlin: 
Springer, 2002) 

P.H. Chavanis, Physica A 387, 5716 (2008) 
Of course, there also exists contracting solutions of the 
Friedmann equations in cosmology corresponding to a 
closed universe [491 ]. 

When d > 2, isothermal distributions and polytropic dis- 
tributions with 7 < 7g/5 = 2d/{d + 2) must be confined 
within a "box" otherwise they have infinite mass. In the 
context of chemotaxis, the box has a clear physical mean- 
ing since it represents the container in which the biolog- 
ical entities live. 

This is not necessary and Eq. (|12|) could be solved in a 
box by decomposing the perturbations on the eigenmodes 
of the Laplacian [Sa ]. 

As is well-known, Einstein [H^ introduced a cosmologi- 
cal constant in the equations of general relativity in or- 
der to obtain a static homogeneous and isotropic uni- 
verse. As a preamble of his paper, he considered a mod- 
ification of the classical Poisson equation in the form 
A"l> — — AirGp. He incorrectly believed that, in 

the Newtonian world model, the cosmological constant is 
equivalent to a screening length. The fact that the cosmo- 
logical constant leads to a Poisson equation of the form 
A"I> = AnGp — A was understood later by Lemaitre and 
Eddington (see J)6, 60, 61] for historical details). On the 
other hand, it was realized that the Einstein universe is 
strongly unstable [g^, so that the idea of a dynamical 
(expanding) universe finally emerged. 
It turns out that the equations of Newtonian cosmol- 
ogy are identical with the Friedmann equations derived 
from the theory of general relativity, provided the pres- 
sure is negligible in comparison with the energy density 
pc^ where c is the speed of light. It is surprising to real- 
ize that Newtonian cosmology was developed after, and 
was influenced by, the theory of general relativity while 
it could have been introduced much earlier. 
The index 74/3 = 2(d— l)/d only appears when we intro- 



but it 



duce the Jeans mass Mj oc p\j oc p('i7-2(d-i))/2^ 
does not play any particular role in the stability analysis 
of the homogeneous fluid. 

In Kostin's approach, the Gross-Pitaevskii equation also 
contains a random potential. In our approach, this 
stochastic term is neglected since we assumed T = 0. 



